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CHAPTER  1 


INTRODUCTION 


1.1  General 


Current  1>  within  the  U.S.  Air  Force  inventory  there  is  cnocgl 
pavement  to  cover  a  two  hundred  foot  vide  section  between  Seattle, 
Washington  and  Miami,  Florida.  A  vast  majority  of  these  pavements  have 
already  reached  the  end  of  their  originally  designed  life.  It  is 
infeasible  for  the  Air  Force  to  consider  replacing  the  aged  pavements 
with  entirely  new  sections  because  the  cost  of  replacing  these  pavements 
to  the  American  taxpayer  would  be  astronomical.  For  this  reason  i;  has 
become  imperative  for  the  Air  Force  to  extend  the  useful  life  of  cost  of 
its  pavements  well  beyond  their  original  design.  To  accomplish  this, 
the  Air  Force  pavement  management  objectives  have  been  to  strengthen 
pavement  sections  where  needed  and  provide  preventative  maintenance  to 
all  sections  in  order  to  forestall  their  evcntwO  functional  failure. 


Kaintcnai’^e  and  rehabilitation  car.  provide  a  good  quality  pavenent 
surface,  is  cost  effective  and  extends  the  serviceable  life  of  the 
pavenent  (1 ,2,3,4,?;. 

The  Air  Force  has  long  recognized  the  need  for  effective  pavenent 
uanagenent  (planning  the  use  and  care  of  pavenents).  In  order  to 
develop  tinely  and  effective  planning,  tne  Air  Force  has  had  to  explore 
rev  and/or  better  netnods  of  creating  an  effective  pavenent  nanagenent 
pin..  This  exploration  was  initiated  in  1970  when  the  Air  Force 
p3\ener.t  responsioilit \  was  assigned  to  a  newly  forced  organization  now 
*  nown  a*-  tu«  Air  Force  Engineer irg  and  Services  Center  (AFESG).  Since 
then  AFfcSC  has  sponsored  and  proooted  research  to  develop  and 
continually  update  an  airfield  pavenent  ranagenent  systec.  The  research 
described  in  this  thesis  is  a  further  devclopaent  of  the  systeas  and 
aids  already  in  use- 

1 .2  Air  Force  Pavenent  Systen  Description 

The  Air  Force  has  hundreds  of  airfields  scattered  through  the 
United  States  and  the  world.  Generally  the  airfields  are  grouped  into 
operational  categories  called  najor  air  cocmnds  (Strategic  Air  Cocxnand 
{SAC},  Military  Airlift  Coanand  (MAC 1 ,  Tactical  Air  Ccocand  ITAC}, 
United  States  Air  Forces  in  Europe  {USAFE},  etc.).  For  the  nost  part, 
these  coanand  groupings  are  not  geographical  but  user  defined. 
Therefore  a  single  aajor  air  coccand  (MAJCOfrJ)  night  *»rcocpass  airfields 
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from  around  the  world  and  at  the  saoe  tine  be  in  close  proximity  with 
another  HAJCOM's  airfields. 

The  decision  structure  fron  vhicn  pavenent  management  planning  is 
developed  is  an  integral  part  of  the  base,  command,  and  Headquarters  Air 
Force  organization.  While  it  is  the  responsbility  of  each  individual 
airfield  (base)  to  develop  a  pavenent  nanagement  plan,  it  is  still 
responsible  to  its  HAJCOM  and  ultimately  the  Air  Force  Headquarters  for 
authority  and  resources  to  execute  it. 


1.3  Long  Range  Planning  in  Pavement  Management 


In  naintaining  an  airfield,  base  pavenent  engineers  must  plan, 
progran,  and  budget  for  maintenance  or  rehabilitative  activities  many 
years  in  advance.  These  engineers  also  seek  assurance  that  their  plans 
provide  the  "best"  available  airfield  for  the  expended  money -  Not  only 
the  needs  of  the  using  aircraft  but  also  the  series  of  rehabilitative 
activities  as  a  whole  must  be  taken  into  consideration.  Only  then  can 
the  pavement  engineers  effectively  manage  the  Air  Force  pavement  system. 

Static,  "What  do  I  do  now?”  approaches  to  pavenent  management  lead 
to  neglect,  waste  and  abuse  of  the  airfield  pavement  system. 
Consideration  of  the  future  as  well  as  existing  status  of  all  airfield 
pavement  sections  and  projected  user  needs  must  be  undertaken. 
Techniques  for  developing  long  range  management  plans  must  provide 
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alternatives  for  decision  makers  to  choose  from,  and  must  be 
sufficiently  flexible  to  allow  input  from  these  decision  makers,  to 
adequately  convey  the  consequences  for  making  a  decision,  and  should 
ideally  direct  the  plan  development  towards  a  goal  or  set  of  goals 
established  by  the  engineers  thenselves. 


1 .4  Research  Objectives 


The  pavenent  nanagenent  needs  of  the  Air  Force  pavement  engineers 
leu  to  the  effort  contracted  to  the  Corps  of  Engineers'  Construction 
Engineering  Research  Lauoratory  (CERL)  in  the  early  1970 's.  As  a 
result,  CERL  has  produced  several  pavenent  nanagenent  aids  (Pavenent 
Condition  Index  [PCI],  Airfield  Pavenent  Management  System  [APMS],  and 
PAVER)  (7,8,97*  These  systems  provide  Air  Force  pavement  engineers  with 
data  storage  and  retrieval,  data  manipulation  and  data  presentation 
capabilities. 

Yet  these  systems  still  constrain  the  pavement  engineers  to  only 
working  in  the  present*  The  designing  and  comparing  of  pavenent 
maintenance  and  rehabilitation  alternatives  remain  directed  at  the 
present  condition  of  the  pavement*  In  order  to  determine  a  "best"  plan 
for  pavement**  for  the  next  twenty  years,  the  engineers  still  have  to 
rely  on  eng-neering  judgement  and  experience.  The  current  systems  do 
not  help  the  engineers  select  and  schedule  pavement  related  procedures 
(overlays,  reconstruction,  etc.)  over  a  period  of  time  (say  20  years). 


These  existing  techniques  do  not  incorporate  optimization  procedures. 
Furthermore  there  is  not  currently  a  way  to  objectively  compare  one 
comprehensive  plan  versus  another. 

Hot  only  are  the  long  range  planning  capabilities  limited  to 
individual  decision  making  for  each  pavement  section,  coordination  is 
lacking  to  optimize  individual  long  range  plans  at  the  network  level  (ar. 
entire  airfield  or  group  of  airfields).  Pavement  engineers  can  not  now 
optimize  the  expenditure  of  funds  for  a  network  of  pavement  sections 
over  a  specified  period  of  time. 

The  research  described  in  this  thesis  has  the  fcl loving  goals: 

1)  The  development  of  a  nethododology  to  optimize  the  selection 
and  timing  of  major  rehabilitative  activities  over  a  specified 
period  of  time  at  a  given  funding  level  for  individual  pavement 
features  (project  level). 

2)  The  development  of  a  nethododology  to  optimize  the  selection  of 
these  activities  at  the  network  level,  also  with  limited 
funding  and  for  a  specified  time  period. 

Simply  put,  develop  a  methodology  to  enable  the  Air  Force  keep 
pavements  in  the  best  condition  possiole  given  a  limited  budget.  As  an 
additional  development,  the  methodology  enables  the  prediction  of  the 
performance  of  the  long  range  plan  as  it  affects  the  whole  network. 


The  research  described  herein  consisted  of  developing  a  methodology 
to  optimize  the  development  of  long  range  airfield  pavement 
rehabilitative  plans,  writing  interactive  user-oriented  computer 
programs  for  implementing  the  methodology,  and  an  analysis  of  the 
information  generated  by  the  programs.  Several  operations  research 
techniques  (dynamic  programming,  integer  programming,  etc.)  were  used 
10  model  the  decision  process  associated  with  long  range  planning  and 
resource  allocation.  An  application  of  the  developer  computer  codes  and 
the  associated  analysis  is  presented  using  an  existing  Air  Force  base. 
The  example  is  explained  in  detail  fr«s  raw  data  gathering  through  an 
discussion  of  the  results. 

Chapter  2  reviews  the  current  Air  Force  pavement  management  systems 
and  discusses  the  ideas,  problems,  and  concepts  associated  with  using 
optimization  in  pavement  management.  Long  range  project  level 
optimization  and  its  related  dynamic  programming  models  are  in  discussed 
Chapters  3  and  4,  while  network  optimization  and  its  integer  programming 
models  are  covered  in  Chapter  3.  An  application  of  the  techniques  and 
analysis  is  made  for  an  Air  Force  base  in  Chapter  6.  Finally  in  Chapter 
7,  the  thesis  is  summarized,  recocendations  are  cade,  and  conclusions 
arc  presented. 
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CHAPTER  2 

PAVEMENT  MANAGEMENT 


2.1  What  is  Pavement  Management? 


Pavement  management  is  a  decision  process  supporting  the 
construction,  maintenance,  and  use  of  aircraft  and  motor  vehicle 
pavements.  motor  vehicles*  According  to  the  Federal  Highway 
Administration  (FHHA),  pavement  management  is: 

"the  process  of  coordinating  and  controlling  all  activities 
related  to  pavements  in  an  attempt  to  best  utilize  public  funds  for 
providing  and  maintaining  pavements  in  a  serviceable  condition  on  a 
continuing  basis.  Effective  pavement  management,  by  necessity, 
involves  the  utilization  of  feedback  of  information  on  pavement 
performance,  pavement  maintenance,  pavement  rehabilitation 

activities,  and  the  cost  of  providing  and  maintaining  pavements" 

(10). 


The  frame  work  or  mechanism  for  pavement  management  is  called  the 
pavement  management  system  (PMS). 

"A  PHS  is  a  set  of  tools  or  methods  that  assist 
decision-makers  in  finding  optimum  strategies  for  providing  and 
maintaining  pavements  in  a  serviceable  condition  over  a  given 
period  of  time.  The  function  of  a  PMS  is  to  improve  the  efficiency 
of  decision-making,  expand  its  scope,  provide  feedback  on  the 
consequences  of  decisions,  facilitate  the  coordination  of 
activities  within  the  agency,  and  ensure  the  consistency  of 
decisions  cade  at  different  management  levels  within  the  sane 
organization"  (11). 

In  general  pavement  management  activities  can  be  characterized  at 
two  levels,  network  and  project.  Network  level  pavement  management 
enables  the  selection  of  the  best  set  of  projects  from  many  sections  of 
the  pavement  system  (i.e.,  state  highway  system,  airfield  or  group  of 
airfields).  At  the  project  level,  decisions  center  upon  specific 
sections  of  pavement. 


The  word  optimization  is  rapidly  becoming  widely  used  in  pavement 
management.  A  clear  understanding  of  optimization,  as  it  applies  to 
pavement  management,  is  essential  before  its  potential  can  be  fully 
realizeo.  The  discipline  of  Operations  Research  (OR)  is  oriented 
towards  finding  beet  solutions  to  problens.  This  search  for  best 
solutions  doesn't  merely  mean  improving  the  status  quo,  but  finding  the 
best  possible  solutions  to  problems.  Usually,  in  order  to  obtain 
solutions,  operations  research  entails  mathematical  modeling  of  the 
decision  process.  But  before  modeling  can  be  undertaken,  there  must  be 
a  clear  understanding  of  the  variables  and  constraints  entering  into  the 
decision  process.  Additionally,  and  equally  important,  the  goals  or 
objectives  of  the  decision  process  must  also  be  understood.  Only  on  the 
basis  of  a  clear  understanding  of  the  relevant  variables,  constraints, 
and  objectives,  can  an  adequate  mathematical  model  of  the  decision 
process  be  constructed. 

A  feasible  solution  to  a  problem  is  one  that  satisfies  all  the 
constraints  in  the  problem.  An  optimal  solution  is  a  feasible  solution 
that  has  the  best  value  of  an  objective  function  (20).  An  objective 
function  is  a  mathematical  representation  of  the  criteria  for  judging 
alternative  feasible  solutions.  Examples  of  objective  functions 
(sometimes  refered  to  as  "benefits")  include  minimize  expended  money, 
maximize  performance,  minimize  the  number  of  people  used,  maximize  total 


traffic  volume,  etc.  In  any  case,  objective  functions  are  used  to 
measure  how  well  a  specific  feasible  solution  satisfied  the  constraints. 
Only  after  all  feasible  solutions  are  explicitly  or  implicitly  evaluated 
can  the  optimal  (best)  solution  be  identified. 

Since  it  is  not  always  possible  to  define,  let  alone  understand  all 
the  variables  that  enter  into  a  problem,  the  optimal  solution  to  a 
mathematical  model  should  be  qualified  as  representing  the  best  solution 
when  considering  only  those  variables,  constraints,  objectives,  and 
structure  of  the  defined  model. 


2.3  Optimization  in  Pavement  Management 


Developments  over  the  last  decade  have  enabled  quantification  of 
the  major  variables  associated  with  pavement  management  (distress, 
performance,  etc).  Due  to  inflated  construction  costs  and  reduced 
pavement  maintenance  budgets,  developing  optimal  decision  policies  has 
been  emphasized.  Although  the  complexity  of  the  pavement  manageaent 
process  precludes  manual  solutions,  the  introduction  of  OR  techniques 
and  computers  has  enabled  effective  analysis.  It  is  now  possible  for 
engineers  to  develop  not  only  feasible,  but  also  optimal  solutions  to 
pavement  management  problems. 

Since  pavement  engineers  seek  to  provide  the  best  possible  pavement 
within  tba  framework  of  many  constraints,  pavement  management  is  a  prime 
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area  for  using  optimisation  techniques.  The  important  attributes  of  the 
pavement  system  (cost,  condition,  performance,  effectiveness, 
etc.)  car.  be  represented  by  quantitative  methods.  Using  these  models, 
pavement  engineers  can  develop  optimal  plans  at  both  the  project  and 
network  levels  of  pavement  management. 

2.3.1  Project  Level  Pavement  Management 


This  level  of  pavement  management  entails  developing  rehabilatative 
strategies  for  specific  sections  of  pavement.  Given  a  feasible  set  of 
alternative  activities  (routine  maintenance,  reconstruction  and  overlay) 
these  activities  should  be  scheduled  in  a  way  that  is  best  with  respect 
to  both  cost  and  benefit  associated  with  the  pavement  section.  Using  OR 
techniques,  alternative  activity  schedules  can  be  developed  for  various 
levels  of  funding.  For  a  specific  level  of  funding,  there  is  at  least 
one  schedule  of  the  alternative  actions  that  provides  the  maximum 
benefit  to  the  pavement  section.  OR  techniques  can  identify  the  best 
alternatives . 


2.3.2  Network  Level  Pavement  Management 


In  addition  to  project  level  pavement  management,  pavement 
engineers  must  develop  a  plan  for  allocating  funds  for  a  group  of 
pavement  sections  (herein  refered  to  as  a  network).  The  network  may 
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represent  only  a  few  sections  (less  than  10)  or  as  many  as  several 
hundred  sections,  and  these  sections  may  be  adjacent  or  discontinous. 
In  any  case,  the  set  of  projects  for  all  pavement  sections  that  provides 
the  most  benefit  oust  be  selected.  For  example,  a  network  analysis  with 
a  hundred  sections,  each  with  three  alternatives,  would  be  difficult  to 
manually  identify  the  group  of  altematves  that  give  the  maximum 
benefit.  But  with  the  aid  of  OR  techniques,  specific  section 
alternatives  can  be  selected  in  order  to  maximize  the  benefit. 


2.3.3  Optimization  Application  in  Long  Range  Planning 


Pavement  engineers  are  faced  with  planning  not  only  for  the  present 
but  also  for  the  future.  If  only  the  present  condition  and  usage  of 
pavements  were  considered  and  if  future  condition  and  usage  is  ignored, 
it  is  impossible  to  maximize  the  structural  performance  of  the  pavement. 
In  developing  a  schedule  of  rehabilitative  measures  over  the  entire 
analysis  period  (say  20  years),  the  consequence  of  scheduling  an 
alternative  in  one  period  significantly  affects  the  scheduling  and 
selection  of  the  alternatives  in  later  periods.  Additionally,  each 
feasible  combination  of  alternative  rehabilitation  techniques  has 
consequences  on  the  performance  of  the  section.  Therefore,  when 
selecting  and  scheduling  a  series  of  rehabilitative  measures,  it  is 
important  to  optimize  the  performance  of  a  pavement  section  over  the 
entire  analysis  period.  Since  the  availability  of  funds  to  support 
pavement  activities  is  highly  variable,  it  is  also  important  to  carry 


out  the  optimization  assuming  several  different  levels  of  funding. 
Because  the  selected  project  schedules  include  work  for  not  only  the 
present  but  also  for  the  future,  choosing  optimally  among  these  project 
alternatives,  at  the  network  level  provides  overall  optimal  long  range 
planning. 


i  Approach  for  this 


The  optimization  approach  in  this  reseach  is  to  initially  develop 
long  range  major  rehabilitative  plans  a^t  the  project  level  for  various 
amounts  of  funding.  These  optimal  plans  for  all  the  pavement  sections 
are  then  used  to  develop  overall  long  range  major  rehabilitative  plans 
at  the  network  level.  Finally,  the  impact  of  various  network  funding 
levels  is  analyzed  to  develop  a  basis  for  estimating  an  appropriate 
level  of  network  funding.  The  methodologies  for  accomplishing  these 
tasks  are  presented  in  the  following  chapters  of  this  thesis. 


2.5  Current  Pavement  1 


.  Systems 


Agencies  responsible  for  pavement  management  whether  local,  state, 
or  federal  utilize  some  form  of  pavement  management  for  planning. 


allocating  budgets,  assigning  work  schedules,  etc.  When  the  pavement 
system  is  not  extensive,  these  functions  can  be  handled  by  a  few 
individuals.  But  as  the  size  of  the  systems  grows,  so  do  the  management 


responsibilities,  problems,  a-*o  stafr  *i/,e  (11). 

The  development  of  pavement  nana^rnent  techniques  has  progressed 
from  the  very  simple  to  the  extrecely  complex,  but  for  tie  purposes  of 
discussion,  four  levels  of  sophisticalioi  will  be  addressed.  In  all 
cases  the  decisions  to  be  made  arc  rich  mairtenance  or  rehabilitative 
activities  are  necessary  end  when  should  they  be  schelduled. 

The  first  level  ot  pavement  management  is  in  essence  r.o  management 
at  all,  that  is  choosing  activities  at  random.  Pavement  engineers  could 
conceivably  cistribute  maintenance  resources  without  any  k**i  fledge  of 
toe  needs  of  either  the  pavement  or  user-  This  is  a  ver\  poor  pavement 
management  system  and  fortunately  rarel>  occures. 

Tne  next  level  ot  sophistication  occures  wren  pavement  engineers 
inspect  the  pavement  network  wit hi:  cheir  jurisdiction  and  establish 
priorities  and  neeos  for  specific  sections  of  the  network.  Priorities 
u  user  requirements  from  either  traffic  composition  and/cr  volumes  are 
specified;  needs  are  determined  from  visual  distress,  technical 
knowledge,  and  the  previous  pavement  experience  of  the  engineers. 
Because  many  pavement  networks  are  complex,  objective  repeatable 
assessments  of  need  are  difficult-  The  subjective  nature  of  these 


decisions  do  not  allow  valid  comparisons  of  alternative  plans. 

The  third  level  of  pavement  management  decision  processes  3Dd 
assessment  systems  are  developed  to  quantify  the  priorities  and  needs  of 
a  pavement  structure.  Objective  management  results,  and  alternatives 
and  user  needs  ere  rationally  evaluated.  At  this  level,  decisions  are 
made  at  the  project  level  and  the  best  activities  selected  considering 
only  individual  sections.  Making  the  mo6t  cost  effective  decisions  for 
each  pavement  section  for  the  benefit  the  entire  network  is  the 
objective  of  the  third  level. 

Today,  the  Air  Force  is  primarily  at  this  third  level  of 
sophistication.  Sophisticated  aids  are  available  (computer  programs 
such  as  APMS,  PAVER,  etc.)  that  allow  pavement  engineers  to 
systematically  develop  objective  alternatives  for  pavements  (8,9).  The 
use  of  computers  has  allowed  examination  of  large  quantities  of  data  and 
detailed  analyses  of  alternative  designs  for  specific  pavement  sections. 
These  programs  use  complex  decision  trees  to  establish  user  and  pavement 
needs,  and  develop  network  plans. 

The  research  in  this  thesis  addresses  a  fourth  level  of  pavement 
management,  the  use  of  OR  techniques  to  develop  optical  pavement 
management  plans.  At  this  level  of  pavement  management,  the  models  and 
procedures  of  the  third  level  are  assumed  to  be  available  and  are 
incorporated  into  an  OR  analysis. 


Without  the  use  of  OR  techniques,  pavement  engineers  have  little 
chance  to  reach  the  "best"  possible  decisions.  Application  of  OR 
techniques  has  been  limited  in  current  pavement  management  systems 
because  of  limited  knowledge  and  only  recent  developments  in  level  three 
management  systems.  The  research  described  in  this  thesis  uses  the 
existing  quantitative  models  of  pavement  performance  and  applies  OR  to 
develop  decisions  for  both  current  and  future  maintenance  and/or 
rehabilitation  activities. 

2 .6  Operations  Research  Applications  in  Paveoent  Management 

Tue  agencies  responsible  for  our  highways  and  airfields  are  at 
various  stages  of  incorporating  pavement  management.  OR  techniques  are 
used  ir  only  a  few  of  the  more  sophisticated  systems,  some  of  which  will 
be  addressed  next. 

2.6.1  Kentucky 


The  Kentucky  Department  of  Transportation  has  developed  methods  of 
applying  dynamic  programming  techniques  to  resurfacing  decisions  to 
optimize  the  allocation  of  expenditures  over  hundreds  of  candidate 
projects  each  year.  A  multistage  d* vision  process  was  proposed  and  the 
optimal  set  of  projects  (based  on  considering  many  pavement  sections, 
each  with  several  alternatives)  could  be  selected. 
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Using  benefit/cost  ratios  the  results  of  dynamic  programming 
generated  solutions  were  compared  to  the  results  of  current  selection 
process.  The  benefit/cost  ratio  achieved  by  the  current  process  was 
only  3.21  compared  to  the  4.22  achieved  by  dynamic  programming 
methodolgy  (19).  The  techniques  have  been  demonstrated  by  the  research 
section  of  the  Department  of  Transportation  but  has  yet  to  be 
implemented  by  the  department  for  routine  use. 


2.6.2  Texas 


The  Texas  Transportation  Institute  has  developed  a  Rehabilitation 
And  Maintenance  System  (RAMS)  for  the  Texas  Department  of  Highways  and 
Public  Transportation.  RAMS  consists  of  seven  computer  programs,  each 
with  its  own  specific  task.  The  combined  and  sequential  use  of  these 
programs  is  expected  to  facilitate  planning,  cost  estimation,  and  fund 
allocation  at  all  levels  of  management  in  the  Texas  Department  of 
Highways  and  Public  Transportation. 

The  OR  techniques  used  in  these  seven  programs  include:  0-1  integer 
linear  programming,  0-1  integer  non-linear  programming,  and  dynamic 
programming.  ‘The  general  objective  of  the  RAMS  programs  is  to  maximize 
the  total  effectiveness  of  all  rehabilitat'  n  and  maintenance  activities 
scheduled  for  the  entire  highway  network  in  the  state  of  Texas  in  each 
year  of  a  predetermined  planning  period  while  remaining  within  the 
available  budget11  (16).  Their  programs  attospt  to  include  not  only 


major  rehabilitation,  but  also  many  other  aspects  of  routine  maintenance 
in  the  decision  process.  Limitations  of  manhours,  equipment,  and 
material  resources  for  the  various  highway  districts  in  the  state  are 
included  as  part  of  the  decision  process.  The  system  has  been 
demonstrated  and  documented  but  has  yet  to  be  implemented  on  a  state 
wide  basis.  Attempts  are  being  made  to  implement  the  system  in  one 
district  next  year. 

2.6.3  Arizona 

The  Arizona  Department  of  Transportation  (ADOT)  has  over  the  past 
several  years  developed  a  network  optimization  system  (KOS)  (by  contract 
to  Woodvard-Clyde  Consultants)  (21).  KOS  uses  linear  programming  to 
solve  for  decision  variables  defined  as  the  proportion  of  roads  in  the 
network  at  the  beginning  of  a  specified  time  period  at  a  specified 
condition  that  will  receive  a  specified  rehabilitation  action.  The 
total  unit  cost  of  these  proportions  is  minimized,  constrained  by 
performance  standards  and  mathematical  limitations  (i.e.,  sum  of 
proportions  be  equal  to  one). 


”It  should  be  noted  that  the  basic  output  ot  the  optimization 
model  does  not  identify  the  rehabilitation  action  that  should  be 
applied  to  a  specific  pavement  in  the  state.  However,  it  is 
possible  to  do  this  by  combining  the  model  output  with  the  data 
base  system  which  stores  road  inventory  data.  Thus,  the  condition 
states  of  various  pavements  in  each  road  category  can  be  identified 
from  the  data  base  system  and  the  results  of  the  optimization  model 
can  then  be  used  to  determine  the  cost  cost-effective 
rehabilitation  action  to  be  applied  to  each  pavement  in  the 
network”  (21). 

The  Arizona  system  is  currently  being  used  on  a  fully  implemented 
basis.  The  systaa  is  still  being  evaluated  for  its  usefulness  and 
validity. 

2.6.4  Ontario,  Canada 

The  Ontario,  Canada  Ministry  of  Transportation  and  Communications 
has  developed  a  pavement  management  system  known  as:  Program  Analysis  of 
Rehabilitation  System  (PARS).  The  computer  based  model  for  planning 
pavement  rehabilitation  programs  for  highway  networks  is  divided  into 
two  main  components:  financial  planning  and  priority  planning. 
Financial  planning  determines  levels  of  funding  for  selected  pavement 
performance  and  usage  goals.  The  second  component,  priority  planning. 
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selects  projects  within  a  fixed  level  of  funding. 

i 

The  selection  and  scheduling  of  projects  for  individual  pavement 
sections  is  done  at  the  network  level  of  analysis.  Linear  programming 
^  models  are  used  to  optimize  the  selection  of  projects  from  a  large  set 

of  feasible  activities.  The  problem  is  formulated  into  a  0-1  integer 
linear  programming  model.  Since  the  problem  is  so  large  (cany  decision 
^  variables)  they  have  relaxed  the  0-1  integer  constraint  and  constrained 

the  decision  variables  to  be  non-negative.  This  allows  the  use  of 
normal  linear  programming  techniques  to  achieve  a  near  optical  solution 
to  the  0-1  integer  model.  The  PARS  system  is  ready  for  use  as  a 
standard  procedure  in  the  Ministry.  The  system  has  not  yet  been 
extensively  used,  therefore  evaluations  of  it  are  not  yet  available. 

I 


2.7  Summary 

This  chapter  discussed  the  process  of  pavement  management,  the 
detimtion  of  optimization  as  used  in  Operations  Research,  and  the 
application  of  optimization  to  pavement  management.  The  next  chapter 
discusses  the  elements  of  an  airfield,  the  project  level  decision 
process,  and  the  adaptation  of  dynamic  programming  to  this  decision 
process. 


CHAPTER  3 


LONG  BADGE  PROJECT  LEVEL  OPTIMIZATION 


3.1  Gg.mil 


This  chapter  discusses  the  decision  process  used  to  develop  long 
range  major  rehabilitative  projects  for  individual  pavenent  sections. 
The  discussion  includes:  explanation  of  airfield  features  (pavenent 
sections),  decision  process  description,  adaption  to  dynamic 
programming,  optimization  criteria,  and  optimization  constraints. 


3.2  Airf ield  Pavenent  Svs^«n  Features 


An  airfield  pavement  system  cannot  be  analyzed  as  a  single  unit 
because  of  the  variability  present  in  pavement  type,  use,  thickness, 
construction  and  maintenance  history,  traffic  ares,  and  condition. 
Therefore  the  airfield  must  be  divided  into  units  with  common 


characteristics  called  "features"  (27).  The  important  characteristics 
of  each  feature  include: 


1) 

Paveaent  type 

2) 

Pavenent  use 

3) 

Paveaent  thickness 

4) 

Construction  history 

5) 

Traffic  area 

6) 

Pavement  condition 

3.2.1  Pavenent  Type 


Paveaent  types  consist  of  flexible,  jointed  concrete,  rigid  overlay 
on  rigid,  nonrigid  overlay  on  rigid,  rigid  on  flexible,  comp  site,  and 
reinforced  rigid  pavements.  A  specific  feature  trust  include  only  one  of 
the  above  types. 


3.2.2  Paveaent  Use 


Airfield  pavenents  are  primarily  used  as  runways,  taxivays,  or 
aprons.  A  feature  must  have  only  one  of  these  uses.  A  taxivay  through 
an  apron  should  be  classified  as  a  different  feature  than  the  apron 
through  which  it  passes. 
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3.2,3  Pavement  Thickness 


Pavement  thickness  varies  greatly  throughout  the  pavement  system. 
Within  each  feature,  the  thickness  should  be  a  constant  nominal 
thickness. 

3.2.4  Construction  History 

On  most  airfields,  the  construction  of  various  portions  of  the 
airfield  occured  at  different  times,  by  different  contractors,  with 
different  materials  and  techniques.  The  construction  and  maintenance 
history  of  the  pavement  within  a  single  feature  must  be  consistent. 


3.2.5  Traffic  Areas 


Airfield  pavements  are  divided  into  “traffic  areas",  based  on  the 
lateral  distribution  of  aircraft  traffic  and  effective  gross  aircraft 
load.  Designated  as  "A",  "B",  "C",  and  "D",  these  areas  represent 
widening  traffic  distribution  and  lessening  gross  loads  (28).  A  feature 
must  be  located  within  a  single  traffic  area. 


3.2.6  Pavement  Condition 

Often  in  an  feature  delineated  with  the  above  criteria,  a  distinct 
portion  of  it  is  in  a  significantly  different  condition.  In  this  case, 
the  area  should  be  redefined  as  a  separate  feature  based  on  the 
differing  pavement  conditions  of  the  two  areas. 

3.2.7  Typical  Airfield  Layout 

Figure  3.1  is  an  example  of  a  typical  airfield  pavement  system.  It 
shows  feature  delineations  depicting  differences  in  pavenent  type,  use, 
thickness,  and  traffic  areas.  Construction  history  is  annotated  as  in 
Table  3.1.  Pavement  condition  is  listed  with  Summary  of  Physical 
Property  Data  table  (see  example  Table  3.2)  under  "General  Condition". 

This  information  is  generally  available  on  location  at  each 
airfield,  in  either  an  airfield  pavement  evaluation  report  or  an 
airfield  pavement  condition  report. 

3-3  Proiect  Level  Definition 

A  long  range  project  level  rehabilitation  plan  represents  the 
rehabilitative  activities  scheduled  on  a  single  airfield  pavement 
feature.  As  alternative  projects  are  developed  for  a  single  feature, 
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ine  operational  interaction  of  th  is  often  not  taken  into  consideration. 
The  project  level  methoct'lcgj  descrioed  in  this  thesis  also  dees  not 
consider  these  interact ic r.s,  because  they  will  be  accounted  fer  at  a 
bigher  level  ot  optimization.  Simply  put,  project  level  optimization  is 
the  development  of  optimal  alternatives  for  a  specific  feature  on  the 
an  f  ield. 

3.4  Pro iect  Level  Decision  Process 

Airfield  pavement  engineers  face  a  complex  decision  process  when 
developing  long  range  project  leve.  major  rehabilitation  plans.  The> 
rust  consider  simultaneously  tnc  i  .ection  and  scheduling  cf  cany 
a mcbil nation  activities,  in  addition  to  the  interaction  of  two  or  more 
activities  and  their  combined  effect  on  the  feature.  The  difficulty  of 
the  problem  is  further  complicated  b,  the  inclusion  of  monetary,  timing, 
and  political  constraints. 

In  this  thesis  the  project  level  oecision  process  is  defined  3S 
selecting  the  best  feasible  activity  at  ^re-established  time  intervals. 
For  each  feature,  the  decisions  are:  Which  rehabilitation  activity 
Mould  be  scheduled  in  each  year  of  the  analysis  period  given  available 
funds.  That  is,  select  the  set  of  activities  that  best  benefit  the 
tavenent  feature  while  at  the  same  time  not  exceeding  a  limited  budget. 
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The  complexity  of  the  decision  process  car.  be  illustrated  by 
considering  the  task  of  scheduling  available  project  activities.  For 
example,  assume  three  activities  (routine  maintenance,  reconstruction, 
and  overlay)  are  available  for  scheduling  at  two  year  intervals  over  the 
next  twenty  years.  In  theory  all  3^®  (59,049)  possible  combinations  of 
these  activities  ‘-hould  be  considered.  If  the  list  of  activities  is 
increased  by  two  possibilites  of  two  thicknesses  of  overlay,  the  number 
of  possible  combinations  to  consider  increases  to  5*®  (9,765,625).  If 
these  activities  must  be  considered  every  year  then  5^  (9.54x10***) 
combinations  must  be  evaluated. 

The  mere  scope  of  this  decision  process  has  previously  prevented 
pavement  engineers  from  considering  all  possibilities.  Figure  3.2 
depicts  the  components  of  this  decision  process. 


3.5 


The  decision  process  previously  described  can  be  modeled  as  a 
staged  decision  process  and  can  be  optically  solved  using  dynamic 
programming,  a  systematic  procedure  for  determining  the  combination  of 
decisions  that  optimizes  the  overall  process.  However,  to  model  and 
solve  the  process,  it  is  important  to  have  a  clear  understanding  of  what 
comprises  the  decision  variables,  the  objective  function,  and  th» 
process  constraints  (291. 
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Decision  variables  are  the  variables  in  the  model  that  can  be 
changea,  i.c.  which  activity  to  choose  for  each  point  in  tirae. 
Changing  the  value  of  a  decision  variable  (selecting  another  activity) 
changes  an  outcome  (performance  of  the  pavement  feature). 

The  objective  function  is  a  measure  of  the  effectiveness  (pavement 
performance)  of  a  specific  set  of  decision  variables  (rehabilitative 
activities).  This  function  is  optimized  i.e.,  the  set  of  decision 
variables  (activities)  is  selected  that  optimize  the  objective  function 
(maximize  the  pavement  performance). 

Constraints  are  mathematical  expressions  representing  conditions 
that  the  decision  variables  must  satisfy.  For  example,  the  total  cost 
of  a  specific  set  of  activities  must  be  within  a  budget  limitation,  the 
condition  of  the  pavement  must  drop  below  a  specified  limit  before  any 
activities  can  be  selected,  etc.  The  decision  variables  must  satisfy 
these  constraints  as  well  as  optimize  the  objective  function. 

With  dynamic  programming,  pavement  engineers  can  solve  the  large 
complex  decision  process  of  selecting  the  best  rehabilitative  activity 
at  specific  time  intervals.  Instead  of  59,049  possible  combinations  to 
analyze  (as  in  the  previous  example),  the  problem  reduces  to  90  possible 
combinations • 
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3.6  Criteria  for  Pavement  Management 


As  previously  noted,  the  objective  function  is  the  mathematical 
expression  that  measures  the  effectiveness  (performance)  of  a  specific 
set  of  decision  variables  (alternative  rehabilitative  activities).  This 
objective  function  becomes  the  optimization  criterion  for  the  decision 
process. 

The  optimization  criterion  used  to  measure  performance  of  a  set  of 
rehabilitative  activities  in  this  thesis  is  the  utility  weighted  area 
under  a  curve  defined  by  the  pavement  condition  index  (PCI)  versus  the 
time  of  the  analysis  period  (7,27). 


3.6.1  Pavement  Condition  Index 


The  pavement  condition  index  (PCI)  is  a  numerical  indicator, 
between  0  and  100,  that  reflects  the  structural  integrity  and  surface 
operational  condition  of  the  pavement-  The  structural  integrity  is  the 
ability  of  the  pavement  to  resist  fracture,  distortion,  and 
disintegration.  The  PCI  is  calculated  as  a  function  of  distress  type 
(rutting,  alligator  cracking,  cracked  slabs,  etc.),  distress  severity 
(depth  of  rut,  spalling  of  pieces  in  alligator  cracking,  number  of  slab 
pieces  and  condition  of  crack,  etc.),  and  distress  density  (extent  which 
distress  covers  sample  area  in  percent  of  total  area).  The  PCI  is  an 
objective  measure  that  closely  correlates  with  the  judgement  of  a  large 
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group  of  experienced  pavement  engineers.  The  method  allows  engineers  or 
technicians  to  rate  neveaent  condition  according  to  a  coo  non  scale  of 
measure. 


To  determine  the  PCI,  the  entire  pavement  systen  (airfield)  is 
first  divided  into  individual  elements  called  features,  as  previously 
defined.  The  -.irfield  is  then  surveyed  feature  by  feature.  The  PCI  can 
be  expressed  mathematically  as  (30): 

p  m. 

PCI  -  100  -  F(t,q)*[  1  Z  a(T  S  =  ,Di:)]  (3.1) 

i=lj=l  3 


where: 


PCI  *=  paveoent  condition  index 
a(  )  *  deduct  weighting  value  depending 

on  distress  type  T.#  level  of  severity 
Sj,  and  density  of  distress 


J 

P 


F(t,q) 


*  counter  for  distress  types 

*  counter  for  severity  levels 

=  total  number  of  distress  types  for 
pavement  type  under  consideration 
=  number  of  severity  levels  on  the 
ith  type  of  distress 

*  an  adjustment  function  for 
multiple  distresses  that  varies 
with  total  suj»ed  deduct  value  (t) 
and  number  of  deducts  (q). 
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3. 6. 1*1  PCI  Det ernination  for  a  Feature  by  Survey 

The  steps  for  performing  the  pavement  survey  and  determining  the 
PCI  for  a  pavement  feature  are  shown  in  Figure  3.3  and  briefly  described 
below  (30): 

1)  The  pavement  feature  is  first  divided  into  sample  units. 

A  sample  unit  for  concrete  pavement  is  approximately  20  slabs; 
a  sample  unit  for  asphalt  is  an  area  of  approximately  5000 
square  feet. 

2)  The  sample  units  are  inspected  and  the  distress  types  and  their 

severity  levels  and  densities  are  recorded.  It  is  imperative 

that  criteria  developed  by  Shahin,  et  al.  (7)  be  used  to 

identify  and  record  the  distress  types. 

3)  For  each  distress  type,  density,  and  severity  level  within  a 
sample  unit,  a  deduct  value  is  determined  from  an  appropriate 
curve  (7).  Step  3  of  Figure  3.3  provides  an  example  of  such  a 
curve. 

4)  The  total  deduct  value  (TDV)  is  determined  by  adding  all  deduct 
values  for  each  distress  condition  observed  for  each  sample 
unit  inspected. 

5)  A  corrected  deduct  value  (CDV)  is  determined  from  the 

appropriate  curve  (7);  the  CDV  is  based  on  the  TDV  and  the 

number  of  distress  conditions  observed  with  individual  deduct 

values  over  five  points  (see  Step  5  of  Figure  3.3). 

6)  The  PCI  for  each  sample  unit  is  calculated  as  follows: 


STEP  6  COMPUTE  PflVEMENT  COMXT10N  NOEX(PCl)lOO*aV  FOR  EACH  SAMPLE  UWT 
KSPECTED 

STEP  7  COMPUTE  PO  OF  ENTIRE  fEATURE  (AVERAGE  PCI'S  OF  SAMPLE  UMTS) 


Figure  3.3  Sumsary  of  Steps  for  PCI  Determination  (27) 
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PCI  =  100  -  CDV  (3.2) 

7)  The  PCI  of  the  entire  feature  is  computed  by  averaging  the  PCIs 
from  all  the  sample  units  inspected. 

8)  The  feature's  overall  condition  rating  is  determined  from 
Figure  3.3,  Step  8,  giving  a  verbal  description  of  the 
pavement's  condition  as  a  function  of  its  PCI  value. 

3. 6 .1.2  Future  PCI  Determination  by  Predictive  Models 

In  addition  to  determining  PCI  of  a  feature  by  direct  survey  of  the 
present  condition,  future  PCIs  of  the  feature  can  be  reasonably 

predicted  with  regression  models.  These  models  were  developed  by  the 
U.S.  Army  Corps  of  Engineers'  Construction  Engineering  Research 
Laboratory  (CERL)  for  predicting  the  consequences  of  different 

maintenance  and  repair  activities  (31) .  Several  initial  models  were 

developed  and  are  currently  being  iuq>roved.  The  models  used  in  this 
research  are  those  from  the  original  work  because  they  are  the  ones 
currently  accepted  and  used  by  the  Air  Force. 

3. 6 .1.2.1  Jointed  Concrete  Pavement  PCI  Prediction  Model 

The  following  model  was  used  to  predict  the  PCI  of  a  jointed 

concrete  surfaced  pavement  feature: 
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PCI  =  100.0  -  AGE* [0. 01967*? AT  -  0.02408*SR 

+  0.001051*(JSL*JSS)  +  2.1057  9* ACOLTHK"0'^ 

+  0 .0347 5*PATCH  +  2.91238  -  0.001775*FI 
-  0.04066*TH4P]  (3.3) 

where: 


?CI  *  Pavement  Condition  Index  at  time  AGE  since 

construction  or  overlay  with  asphalt  or  concrete 
AGE  **  tine  since  construction  of  slab  or,  if  overlaid,  tine 
since  overlay  construction  (years) 

FAT  *  (ratio  of  interior  slab  stress/oodulus 
of  rupture)  x  100 

SR  =  slab  replacement  (percent  total  slabs) 

JSL  *  longest  joint  spacing  (feet) 

JSS  »  shortest  joint  spacing  (feet) 

ACOLTHK  =  thickness  of  overlay  (inches) 

FATCH  =  slabs  containing  large  patches 

(5  square  feet),  percent  of  total  slabs,  or 
percent  area  of  total  area  patched  if  overlaid  with 
asphalt 

TEMP  =  average  annual  teaperature  (°F) 

FI  =  freezing  index  (degree  days  below  32°F) 

Note,  when  the  pavement  is  not  overlaid,  the  entire  ACOLTHK  term  is 
eliminated  because  it  is  undefined  when  it  is  zero. 


The  difference  between  this  Equation  3.3  and  the  original  one 
derived  by  CERL  is  in  the  term  ACOLTHK  and  therefore  does  not  reflect 
the  sane  prediction  characteristics  as  the  CERL  equation.  The  CERL 


equation  had  a  binary  (0  or  1)  ACOLTHK  variable.  Hence  it  was 
insensitive  to  the  various  thicknesses  of  asphalt  overlays;  and  the 
resulting  predicted  PCI  values  were  not  realistic.  The  new  term  was 
derived  fron  a  data  set  of  31  points  and  plotted  as  shown  in  Figure  3.4. 
The  original  equation  was  based  on  an  average  overlay  thickness  of  2.7 
inches.  The  resulting  coefficient  of  the  binary  term  was  0.94191.  The 
new  tern  ACOLTHK  (overlay  thickness)  gives  the  same  result  when  2.7 
inches  of  overlay  are  used,  but  makes  the  total  PCI  predictive  equation 
much  more  sensitive  to  changes  to  the  overlay  thickness.  Figure  3.5 
shows  the  sensitivity  of  the  two  equations  as  a  function  of  time.  Both 
equations  are  plotted  for  2,  4,  and  6  inches  of  overlay.  Note  the 
change  in  spread  of  the  two  equations  as  a  function  of  thickness. 

For  this  equation  to  be  useful  in  the  dynamic  programming  model, 
all  the  terns  must  be  calcuiaDie  whenever  the  PCI  is  needed  in  the 
analysis.  The  methods  used  to  determine  the  values  for  these  variables, 
are  discussed  in  the  next  chapter  under  model  development. 

3. 6 .1.2. 2  Asphaltic  Concrete  Surfaced  Pavement  PCI  Prediction  Model 

The  following  model  was  used  to  predict  the  PCI  of  an  asphalt 
surfaced  pavement  feature  (31) : 

PCI  =  100  -  AGE*ll.487/.isg  +  0.143*AGECOL  +  6.56/Tac  -  1.23*aac]  (3.4) 


Figure  3.4  Change  in  PCI  per  Overlay  Thickncfl# 
Verbus  Ovei lay  Thickness 


100 


Figure  3.5  Example  Showing  Sensitivity  of  PCI  Equating 
to  Overlay  Thickness 


AGE 


age  since  original  const rt-cr ion  or  ='.cce  last 
overlay  if  the  pavement  has  been  cverleid 


a 


sg 


AGECOL 


lead  repetition  factor  determined  at  t  .a  ''abgrace 
level;  sg  is  a  function  of  total  pavement 
thickness  above  the  subgrade,  subgrade  CBR,  and 
the  tire  contact  area  and  tire  pressure  of  an 
equivalent  single  wheel 

age  betveen  the  tine  the  pavement  was  constructed 
and  the  tine  it  received  the  last  overlay;  equals 
zero  if  the  pavem?nt  was  net  overlnd 
total  asphalt  thickness  in  inches  including  overlay 
load  repetition  factor  determined  at  the  asphalt  base 


The  assignment  of  value-  to  the  terms  in  this  equation  is  3lso  discussed 
m  the  next  chapter  under  model  development. 


2,6.2  PCI  Utility 


There  is  a  difference  in  value  to  pavement  engineers  in  improving  a 
pavement  from  different  levels  of  PCI.  For  example,  improving  the  PCI 
of  pavement  from  95  to  100  is  far  less  significant  than  improving 
another  paveuent  from  50  to  55.  This  crangc  in  value  for  a  change  in 
PCI  as  a  function  of  PCI  can  be  modeled  with  utility  theory.  Utility  is 
the  true  measure  of  value  to  pavement  engineers. 


I.- roups  of  txpt  r.encc-d  .Ur  Force  pavement  engineers  were  used  to 
ostab1is.li  the  utility  cl  the  PC  i  as  a  function  of  pavement,  type  (primary 
er  secondary).  The  results  of  this  survey  are  shown  in  Figures  3.6,  3.7 
‘•uf  -5*8  (32).  These  scales  reflect  PCI  utility  as  a  function  of  user 
•eeds  and  ir.  this  rese.aich  usee  needs  art1  accounted  for  in  the  network 
iuv'el.  analysis  discussed  in  later  chapters.  Therefore  a  PCI  utility 
mode  L  war.  extrapolated  from  the  maintenance  and  repair  zone  versus  PCI 
ible  developed  by  a  group  of  experience  Air  Force  pavement  engineers 
(Figure  3.9). 


Using  both,  the  utility  curves  (Figures  3.6,  3.7  and  3.3  and  the 
correlation  of  maintenance  and  repair  (Figure  3.9),  the  utiLity  of  the 
’ll  was  developed  by  the  author  foi  the  purposes  of  this  research 
‘.figure  3.10).  The  breaks  in  the  curve  at  70  and  25  PCI  represent  the 
inclusion  of  major  rehabilitation  activities  for  consideration  and  the 
;  .auditory  use  of  overall  rehabilitation  respectively.  The  linear 
i  id-portion  of  curve  closely  resembles  the  shapes  of  the  previously 
developed  curves  (Figures  3.6,  3.7  and  3.8). 


Integrating  this  curve  and  plotting  the  values  of  the  integral  with 
respect  to  PCI  results  in  Figure  3.11.  This  utility  weighted  PCI 
reflects  the  engineers'  bias  in  relative  value  of  FCI  as  a  function  of 
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PAVEMENT  CONDITION  INDEX 


Figure  3.7  Utility  Curves  for  Taxiways  (32) 
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PAVEMENT  CONDITION  INDEX 

Figure  3.8  Utility  Curves  for  Aprons  (32) 


UTILITY 


>6 


M  &  R  ZONE 

PCI 

!  RATING 

ROUTINE 

IOO 

up 

Us 

85 

VERY  GOOD 

_ i 

ROUTINE, 

MAJOR, 

OVERALL, 

70 

■ 

55 

||  |  pi 

ill  ill  111 

X 

|  FAIR 

MAJOR, 

OVERALL 

40 

n 

| 

29  POOR 

m 

& 

OVERALL 

25 

j-il  VERY  POOR 

10 

0 

& 

FAILED 

Figure  3.9  Correlation  of  fcaintenance  and  Repair  Zones 
and  PCI  and  Condition  Rating  (32) 
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Utility  o£  PCI  Curve 


Figure  3.11  Utility  Weighted  PCI  Vcrsuu  PCI 
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3.6.3  Rehabilitative  Activity  Performance  Measure 


Using  PCI  and  its  utility,  the  performance  of  a  specific 
rehabilitative  activity  can  be  measured  over  time.  The  performance  can 
be  measured  as  the  area  under  a  utility  weighted  PCI  versus  time  plot 
(Figure  3.12).  If  a  rehabilitative  activity  is  applied  to  a  pavement 
(as  indicated  by  the  vertical  spike  in  the  curve  on  Figure  3.13),  the 
performance  of  the  activity  is  measured  by  the  area  under  the  curve.  A 
different  activity  results  in  a  different  area.  The  area  can  be 

calculated  by  integrating  of  the  equation  involving  the  PCI  predictive 
equation  and  the  PCI  utility  function  equation,  as  follows. 

The  utility  weighted  PCI,  as  shown  in  Figure  3.11,  can  be 

represented  by: 

°PCI  =  3  +  +  c*PCI2  +  d*PCI3  (3.5) 

where  a,  b,  c,  and  d  are  the  regression  coefficients: 
a  **  -1.59 
b  *  1.33 
c  *=  -0.01 
d  *  1.17  x  NT5 

The  performance  of  the  pavement  feature  as  defined  by  the  area  under  the 

utility  weighted  PCI  versus  time  plot  can  be  solved  by  integrating  the 

utility  of  PCI  with  respect  to  age  (A): 


°PCI  dA 


Performance  = 


(3.6) 


TIME,  YEARS 


Figure  3.12  Utility  Weighted  PCI  Versus  Tice  Curve 


lOd  031HSI3M  AIHlin 


Figure  3.13  Performance  of  Rehabilitative  Activity 
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A  =  AGE  =  age  of  pavement  feature  in  years. 
By  substituting  Equation  3.5  into  3.6  we  get: 


Performance  =  f  (a  +  b*PCI  +  c*PCI^  +  d*PCI^)  dA  (3.7) 


Depending  on  whether  the  feature  is  a  Portland  concrete  surfaced  or  an 
asphaltic  concrete  surfaced  pavement,  the  appropriate  PCI  predictive 
equation,  3.3  or  3.4,  is  substituted  into  Equation  3.7.  In  general 
Equations  3.3  and  3.4  can  be  represented  as: 


PCI  =  100  -  AGE  *  K 


K  =  parameter  representing  the  structural,  climate, 
and  maintenance  characteristics  of  the  feature, 
either  asphalt  or  concrete  surfaced. 
Substituting  Equation  3.8  into  Equation  3.7, 


Performance  =  ’  la  +  b(100-AGE*K)  +  c(100-AGE*K)^ 

+  d(100-AGE*K)3l  dA 

and  evaluating  Equation  3.9  evaluated  from  time  to  time  T2  as: 
T 

Perfornance  L2=a*AGE  -  b(100-AGE*K)2/2X  -  c(100-AGE*K)3/6K 

1  ,  T2 
-  d(100-ACE*X)4/12K  |  ' 


With  the  appropriate  substitution  for  K,  Equation  3.10  can  be  used  to 
easily  calculate  the  area  under  the  utility  weighted  PCI  versus  tine 
plot  between  time  T ^  and  T2  in  years. 

3.7  Pavement  Management  Constraints 

There  are  constraints  that  must  be  considered  when  developing  long 
range  rehabilitative  plans.  In  order  for  the  model  to  resemble>  as 
closely  a6  possible,  the  engineers'  decision  process,  these  constraints 
must  be  included  in  the  dynamic  programming  model.  The  following 
constraints  were  incorporated  in  the  model  in  this  research: 

1)  Budget  licit  -  total  amount  which  the  plan  must  not 

exceed  (in  present  worth  dollars). 

2)  Tine  in  years  -  total  time  plan  must  not  exceed  (in  years) 

3)  Decision  intervals  -  scheduled  points  in  time  where 
activities  are  considered. 

A)  Minimum  Distress  -  level  of  PCI  which  pavement  must  reach 
before  any  rehabilitation  can  be  considered. 

5)  Limitation  in  Activity  Selection  -  some  alternatives  can  not 

be  repeated  (i.e.,  reconstruction)  during  the  analysis  period. 
Also  overlays  must  be  controlled  by  limiting  the  total  amount 
of  asphaltic  concrete  applied  during  the  analysis  period. 
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3.8  Summary 


This  chapter  discussed  the  elements  of  an  airfield,  the  decision 
process  at  the  project  level,  adaptation  of  dynamic  programming  to  the 
decision  process,  and  the  optimization  criterion  and  constraints 
involved  in  the  process.  The  next  chapter  expands  on  these  topics  as 
the  components  of  the  computer  program  developed  in  this  research,  for 
the  project  level  decision  process  is  explained. 
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CHAPTER  4 

DYNAMIC  PROGRAMMING  MODEL  FOR  PROJECT  LEVEL  OPTIMIZATION 


4.1  General 


This  chapter  discusses  the  application  of  dynamic  programming 
techniques  to  long  range  project  level  optimisation.  The  discussion 
includes:  explanation  of  the  dynamic  programming  method  and  the 
mechanics  of  applying  it  to  the  project  level  decision  process. 


4.2  Dynamic  Programming  Method 


Dynamic  programming  is  an  approach  to  optimization  vhich  is  useful 
in  solving  specific  kinds  of  problems.  It  is  not  a  particular  algorithm 
in  that  it  is  not  a  specific  procedure  or  set  of  rules  for  finding  the 
optimal  set  or  solution. 
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"Dynamic  programming  is  a  way  of  looking  at  a  problem  which 
nay  contain  a  large  nunber  of  interrelated  decision  variables  so 
that  the  problem  is  regarded  as  if  it  consisted  of  a  sequence  of 
problems,  each  of  which  required  the  determination  of  only  one  (or 
a  few)  variables.  Ideally,  what  we  seek  to  do  is,  in  effect, 
substitute  solving  n  single  variable  problems  for  solving  one  n 
variable  problem.  Whenever  this  is  possible,  it  usually  requires 
very  mich  less  computational  effort.  Solving  n  smaller  problems 
requires  a  computational  effort  which  is  proportional  to  n,  the 
number  of  single  variable  problems  if  each  problem  contains  one 
variable.  On  the  other  hand,  solving  one  larger  problea  with  n 
variables  usually  requires  a  computational  effort  which  is  very 
roughly  proportional  to  a11,  where  a  is  soae  constant.  Hence  the 
desirability  of  transforming  or  considering  an  n-dimensional 
problem  as  n  one-dimensional  problems"  (29). 


This  method  was  first  introduced  by  Bellman  (41)  with  a  stated 
principle  ot  optimality.  When  applied  to  the  project  level  decision 
process,  Bellman's  principle  simplifies  to: 


Any  optical  long  range  rehabilitation  plan,  consisting  of 
major  rehabilitative  alternatives  scheduled  at  selected  points  in 
the  analysis  period,  is  comprised  of  a  sequence  of  smaller  optimal 
sub -plans. 


For  example,  the  activities  scheduled  by  an  optimal  20  year 
rehabilitation  plan  prior  to  the  10  year  point  are  still  optical  if  the 
analysis  period  is  reduced  to  first  10  years. 

The  first  advantage  of  dynamic  programming  has  already  been 
mentioned:  the  required  search  in  a  coeplex  decision  tree  frca  a  single 
n-dinensional  optimization  problen  is  reduced  to  n  one-dicensional 
probleos,  which  can  be  solved  one  at  a  time.  A  second  advantage  of 
dynacic  programming  is  that,  within  the  framework  of  the  models  used, 
global  optical  solutions  can  be  reached.  Third,  the  eabedded-solut ion 
nature  of  dynanic  programing  is  useful  because  the  saae  problem  does 
not  need  to  be  repeatedly  solved  for  various  analysis  period  lengths. 
Solving  the  problec  for  the  longest  useful  analysis  period  under 
consideration  provides  results  for  all  analysis  periods  less  than  the 
original.  Hence  considerable  cocputation  tine  nay  be  saved.  Finally, 
although  decision  variables  that  are  required  to  be  integer  or  discrete 
create  difficulties  for  other  optinization  techniques,  their  presence 
actually  simplifies  the  dynanic  programing  solution  process. 


The  project  level  decision  process  can  be  depicted  as  a  complex 
decision  tree  with  a  choice  of  five  activities  (feasible  decisions)  at 
pre-established  periods  of  tine  (stages).  For  example,  Figure  4.1 
illustrates  the  initial  state  of  the  pavenent  feature  at  tine  T^  ancj 
three  succeeding  tine  intervals  T^,  T2,  and  T3.  Suppose  that  five 
possible  activities  are  available  at  each  of  these  tine  intervals. 
After  three  intervals  there  are  5^  (125)  possible  combinations-  If  the 
decision  tree  were  extended  to  ten  intervals,  a  total  of  5*®  (9,765,625) 
possible  paths  would  exist.  If  a  coaputer  took  0.1  second  to  calculate 
the  necessary  data  for  analysis  of  each  path,  it  would  require  nore  than 
113  days  to  enuaerate  all  paths.  Using  dynamic  programing  isetbods,  the 
sane  problen,  over  a  ten  interval  analysis  period,  could  be  solved  in 
less  than  25  coaputer  seconds  and  still  insure  optinality. 

Figure  4.2  is  the  decision  tree  for  the  sane  process  as  in  Figure 
4.1  except  with  dynaaic  programing  nethods  applied.  Note  that  the 
nucber  of  possible  states  (condition  of  pavement)  remains  the  saoe  (5) 
each  tine  the  five  activities  are  applied.  The  five  decisions 
(rehabilitative  activities)  are  applied  to  each  of  the  five  entering 
states  (pavenent  conditions).  The  best  decision  (rehabilitative 
activity)  is  selected  for  each  o£  entering  five  conditions.  These  "best 
activity  selections"  reduce  the  number  of  resulting  states  to  five. 
This  process  is  repeated  for  each  stage  (analysis  time  period).  At  the 
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end  of  the  analysis  period  with  the  activity  selections  oade  for  each  of 
the  entering  five  states,  another  selection  of  these  five  states  is  made 
to  determine  the  best.  Selecting  the  best  of  these  five  terminal 
decision  paths  (states)  also  establishes  the  overall  (global)  best 
(optimal)  path  (set  of  activities).  The  activities  selected  are  found 
by  retracing  the  optimal  path  backwards  through  the  tree. 

Hot  only  does  the  dynamic  programming  method  reduce  the  mraber  of 
paths  to  be  analyzed,  but  it  also  reduces  the  amount  of  computer  memory 
necessary  to  store  the  number  of  paths  which  must  be  saved  for  possible 
retracing. 

Summarizing,  a  decision  path  is  set  of  rehabilitative  activities 
scheduled  at  specific  tine  intervals  during  the  analysis  period. 
Dynamic  programming  methods  reduce  the  total  number  of  decision  paths  to 
be  analyzed,  searching  for  the  optimal  or  best  path.  Traditional 
decision  tree  analysis  is  multiplicative  in  nature  as  the  cumber  of 
stages  (analysis  periods)  increases.  Whereas,  with  dynamic  programming 
methods,  the  decision  process  is  only  additive  in  nature. 


4.3.1  Dynamic  Programming  Model  Flow  Chart 


The  dynamic  programming  method,  as  applied 
decision  process,  is  illustrated  in  Figure  4.3. 


to  the  project  level 
Each  stage  r presents  a 
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decisions  (activities)  applied  at  each  state  are: 


1) 

Routine  Maintenance  (R/M) 

2) 

Reconstruction  with  PCC  (R/C) 

3) 

Overlay  with  2  inches  of  AC 

(0/L2) 

4) 

Overlay  with  4  inches  of  AC 

(0/L4) 

5) 

Overlay  with  6  inches  of  AC 

(0/L6) 

Substitution  or  addition  of  other  activities  could  be  done.  These 
decisions  are  represented  as  the  vector  [D)  and  defined  as  follows: 

ID)  «  Idj,  d2,  d3,  d4,  d5)  (4.1) 

The  single  independent  state  variable,  PCI,  is  required  for  this  dynamic 
programming  formulation.  When  entering  into  any  stage  i,  it  can  have 
only  one  of  five  values.  Associated  with  the  PCI  are  14  other  dependent 
state  variables. 

Each  of  the  five  possible  entering  states  is  defined  by  a  single 
row  vector  of  independent  and  dependent  state  variables  values  as: 

lsij'  “  taU«  ai2’  ai3>  (4*2) 

with  a^  as  one  of  the  following  attributes: 

1)  Pavement  Condition  Index  (PCI) 

2)  Performance  or  area  under  curve  (BENEFIT) 
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3)  Expended  cost  (COST) 

4)  Age  of  the  pavenent  (AGE) 

5)  Alpha  under  the  asphalt  layer  (aflc) 

6)  Alpha  on  the  subgrade  (cigg) 

7)  Age  of  the  overlay  (AGECOL) 

8)  Thickness  of  the  asphalt  layer  (T  ) 

9)  Thickness  of  the  PCC  layer  (T  ) 

J  pcc' 

10)  Percent  slabs  replaced  (SR) 

11)  Whether  or  not  the  pavenent  has  been  overlaid  (ACOL) 

12)  Percent  of  area  patched  (PATCH) 

13)  Age  since  original  construction  (IAGE) 

14)  Whether  flexible  or  rigid  pavenent  (PVTTYP) 

13)  Selected  rehabilitative  activity  (DECIDX) 


i  **  index  of  the  possible  enter  states 
j  =  index  of  the  attributes 


The  entire  entering  set  of  states  for  the  pavenent  can  be 
represented  in  a  single  array,  [Sj  for  each  stage  n  and  defined  as: 


all»  a12*  a13»  •••*  aij 
a21 »  a22>  a23>  •••»  a2j 


31 »  a32»  a33»  •••>  a3j 


ij>  ai2*  ai3»  •••*  3ijJ  n 
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where:  n  =  index  of  stages 

The  return  value  (optimization  criteria)  is  used  to  select  the  best 
activity.  The  return  function  in  this  model  is  the  total  area  under  the 
utility  weighted  PCI  versus  time  curve.  The  values  of  the  return 
function  are  represented  by  [ R ]  and  defined  as: 

lR]  =  Irj,  r2>  r^»  r^,  r^)  (4.4) 

where  rQ  represent  the  return  of  the  selected  activity  for  each  state 
when  all  the  activities  are  applied  to  each  entering  state.  The  optical 
(best)  returns  are  then  kept  and  carried  forward  to  the  next  state  as 
cenbers  if  the  st3te  variable  array,  (Sln.  They  are  stored  as  the 
performance  attribute,  BENEFIT. 

The  dynamic  prograntsing  model  can  be  represented  na theca tical ly  as 
follows : 

Maximize: 

n  T 

If  l-1.59+1.33*PCI._o.01*PCI?+i.l7*10_5*PCl3)  dtine  (4.5) 

L*1  Tx  1 

with  the  following  constraints  or  attribute  transformations: 

PCI  limits  for  all  stages  _i 
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0  £PC1i  <100  for  all  i  (4.6) 

Determination  of  PCI  for  al  1  staees  _i 


PClj^  =  ioo  -  AGE100  *  K 


Decision . 


feasibility  for  all  stages  j. 


(4.7) 


If  0  <  PCI <_  70  then  the  feasible  decisions  are: 
dl>  d2>  d3>  d4»  d5 

Else  if  70  <  PCI.  j  100  then  the  only  feasible  decision  is: 

dj,  routine  maintenance  only  (4.8) 

Project  cost  1 imitation  for  all  stages  i 


COST^  <  Project  Funding  Limit  (MAXCOST)  (4.9) 

Age  of  the  pavement  for  all  stages _i 

AGE^  =  AGE^_j  +  stage  i  interval  length  (4.10) 

Thickness  if.  asphalt  layer  for  all  stages  JL 
If  the  decision  at  the  beginning  of  stage  i  was: 


dl  then  Tac^  =  Tac^_j 
d2  then  Tac^  =  0 
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d3  then  Tac^  =  Tacj.j  +  2 
d4  then  Tac^  =  TaCj_j  +  4 

d5  then  Tac^  =  TaCj_j  +  6  (4.11) 

Thickness  of  Portland  cement  concrete  slab  for  al 1  stages  _i 
If  decision  at  the  beginning  of  stage  i  vas : 

^l»d3>d4»d5  then  Tpcc^  =  Tpcc^j 

then  Tpcc^  =  reconstructed  pcc  thickness  (4.12) 

Structural  parameters  of  flexible  pavement  for  al 1  stages  _i 

ai  =  (Ti/A0-5)il0.048-1.15621og(CBR/pe) 

-0.064141og(CBR/pe)2-0-4731og(CBP./pe)3]  (4.13) 

where : 

=  thickness  of  layer  above  point  of  interest  in  pavement 
(just  under  AC  and  on  top  of  the  subgrade) 

A  *  contact  area  of  one  tire  of  aircraft  (psi) 

CBR  *  California  bearing  Ratio  of  underlying  layer 
pe  =  tire  pressure  (psi)  calculated  using  contact 

area  A,  and  the  equivalent  single-wheel  load  (ESWL) 
determined  at  depth 
Age  ot  overlay  for  al  1  stages  _i 
If  the  decision  at  the  beginning  of  stage  i  was: 


d!  then  AGECOL  ^  »=  AGECOL^j 
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<*2  then  AGECOLj  =  0 

d3,d4,d5  then  AGECOL^  =  IAGE^  (4.14) 

Percent  of  slabs  replaced  for  al 1  stages  _i 

SRi  =  88.622-1. 993*PCIi+O.012*PCI?-1.127*PCI?  (4.15) 

Pavement  overlay  status  for  al 1  stages  _i 

If  the  decision  at  the  beginning  of  stage  i  was: 

dj  then  ACOL^  =  ACOLi_1 
dj  then  ACOL^  =  0 

d3.d4.d5  then  ACOLj  =  1  (4.16) 

Percent  of  feature  patched  for  al 1  stages  i 

PATCH ^  =  32. 204+0. 1306 *PCIi  -0.0098*PClf  +5.28*10-5  *PCI?  (4.17) 

Age  since  original  construction  or  reconstruction  for  al  1  stages  i_ 

If  decision  at  the  beginning  of  stage  i  was: 

^1»^3-*^A»^5  then  lAGE^  =*  IAGEj^j.  +  stage  interval  length 

d2  then  IACEi  =  0  (4.18) 

Pavement  type  index  (rigid  =  1  and  flexible  =  2  )  for  al  1  stages  J. 

If  the  decision  at  the  beginning  of  stage  i  was: 


dl,d3»d4>d5  then  PVTTYP^  =  PVmfPi,! 


&2  then  PVTTYP^  ■=  1 

Selected  rehabilitation  activity  for  al 1  stages  X 
If  the  decision  at  the  beginning  of  stage  i  was: 

<4  then  DECIBXj  =  1 
<*2  then  DFCIK^  =  2 
d3  then  DECIDXj  -  3 
then  DECIDX^  =  4 
d5  then  DECIDX j  =  5 


(4.19) 


(4.20) 


where: 

d  ,d2,...,d5  c  decision  of  routine  maintenance jreconstruct ion 
with  pcc;  and  overlays  of  2,4, 
and  6  inches  of  AC  respectively 
i  *  stage  index 

n  *=  number  of  stages  in  analysis  period 
AGE100  «*  time  since  pavement  was  at  PCI  *  100 
If  the  feature  was  rigid  then: 


K  «  0 . 01 96 7*F AT-0 . 0240 8*SR+0. 001051  *(JSL*JSS) +2. 1057  9*AC0LXfiX“°*d^ 

■K)  .0347  5*PATCH+2 .  91 23  S-0  .00177  5*FI-0 .0406  6*TEMP  (4.21) 


Else  if  the  feature  was  flexible  then: 


K-1.48//  .sg40.143*AGEC0L^.56/Tac-1.73*  o3c  (4.22) 
Note  that  the  variables  needed  to  determine  K  are  explained  in  Chapter  3. 

The  single  independent  state  variable  in  this  dynamic  programming 
algorithm  is  PCI.  Associated  with  this  state  variable  is  a  vector  of 
state  variable  attributes  dependent  on  the  PCI  state  variable  and  the 
feasible  decisions.  Because  of  the  nature  of  the  decisions  available, 
at  any  stage  (decision  point)  during  the  analysis  period  there  can  be  no 
nore  than  five  values  of  the  state  variable. 

4.3.2  Computer  Program  Description 

A  computer  program  was  written  for  using  dynamic  programming 
methods  for  selecting  and  scheduling  the  set  cf  rehabilitative 
activities  that  optimize  (maximize)  the  objective  function  (performance 
as  defined  as  the  utility  weighted  PCI  versus  tine  curve). 

The  program  was  written  in  Control  Data  Corporation's  Fortran  V 
(which  generally  adheres  to  the  Fortran  77  ANSI  standards).  It  is  a 
2900  line  program  with  73  subroutines  supporting  the  ^ain  program. 
Figure  4.4  is  a  flow  chart  showing  the  major  elements  of  tb  .  dynamic 
programming  algorithm  as  applied  to  the  project  level  decision  process. 
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After  the  feature  data  and  the  analysis  criteria  are  entered  via 
condensed  data  secs,  the  prcblen  is  initialized.  As  shown  in  the  flow 
chart,  figure  4.4,  the  decision  process  is  repetitive,  applying  the  five 
decisions  to  each  of  the  five  states  at  eacn  stage,  and  generating  the 
five  new  states  for  the  next  stage.  The  process  is  continued  until  all 
the  stages  in  the  analysis  period  are  ccspleted. 


4.3.3  Coaputer  Progran  Inputs 


Input  data  for  the  computer  progrsa  are  read  frer.  two  condensed 
data  files:  analysis  criteria  data  and  feature  property  data. 


4. 3. 3.1  Analysis  Criteria  Data 


The  data  necessary  for  characterizing  the  dyuanic  progressing 
analysis  are  input  fron  a  separate  data  file.  The  information  are  r tsA 
in  condensed  fora  as  follows: 

1)  Nunber  of  stages  or  intervals  in  the  analysis  period 
(NUKSTG) 

2)  Ages  o t  each  interval  in  analysis  period,  stage  ages  (STGAGE) 
and  terainal  age  of  analysis  period,  final  stage  age  (FSTGAGE). 

3)  Condition  rating  of  pavenont  (PCI)  above  which  ac 

rehabilitative  activities  will  be  accomplished  (PCXHIH) 
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A)  Intere°t  and  inflation  rates  for  use  during  the  analysis  period 
(INTRST  and  IHFLTN) 

3)  Control  variables  to  restrict  the  type  and  amount  of  output. 

A-3.3.2  Pavement  Feature  Data  Input 


A  condensed  data  set  is  read  into  the  program  free  a  separate  data 
file  containing  features  in  the  network.  After  finishing  the  analysis 
on  one  pavement  feature,  tht  program  reinitializes  the  entire  process 
and  reads  in  data  for  another  feature  and  continues  until  all  have  been 
analyzed.  The  data  read  into  the  program  for  each  feature  are: 

1)  Base  or  airfield  name  (BNAHE) 

2)  Pavement  feature  name  (FNAHE) 

3)  Pavement  type,  rigid  or  flexible  (PTYPE) 

4)  Feature  length,  width  and  current  age  (LNGTH,  WDTH  and 
LACE) 

5)  If  the  feature  is  flexible  then,  number  of  layers  in  pavement, 
overlay  thickness,  age  of  the  overlay,  and  the  physical 
properties  of  each  layer,  thickness,  type  and  CBR  (NUHBKLYR, 
OATHK,  AGECOL,  LYBXYP,  LTRTEK,  and  LYRCBR).  Additionally,  the 
current  PCI,  average  annual  tesperature,  and  freezing  index  are 
listed  (OIPCI,  TcMP,  and  FI). 

6)  If  the  feature  is  rigid,  the  pavement  characteristics  are 
input,  slab  size,  modulus  of  rupture,  modulus  of  soil  reaction. 
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slab  thickness,  overlay  thickness,  current  PCI,  percent  slabs 
replaced,  freezing  index,  percent  of  feature  area  patched, 
average  annual  tenprature,  and  the  age  of  the  last  overlay 
(JS1,  JSS,  KDRPTR,  KMPD,  OPTHK,  OATHK,  01  PCI,  SR,  FI,  PATCH, 
TEMP,  AGECOL). 

7)  Design  aircraft  data,  identification  and  weight  (ARCFTIB  and 
ARCFTWT) . 

8)  Using  aircraft  type,  mission,  and  number  (ARCFTCD,  ARCFTKS,  and 
ASCFTH ) . 

9)  Pavement  feature  use  2nd  need  (PVTUSE  and  PVTND). 

10)  Feasible  rehabilitative  activity  data,  maximum  asphaltic 
concrete  overlay  thickness,  reconstruction  Portland  cement 
concrete  thickness,  overlay  cost,  reconstruction  cost, 
reconstructed  eodu)us  of  rupture,  reconstructed  modulus  of  soil 
reaction,  reconstructed  slab  size  (HXACTBK,  RPCISK,  GLCST, 
RCSTCST,  RHDRPTR,  RKHOD,  RJSL,  2JSS) 

4.3.4  Computer  Progisa  Algorithms 


In  addition  to  the  algorithms  used  to  describe  the  PCI,  PCI 
utility,  and  the  area  under  the  utility  weighted  PCI  versus  time  curve, 
functions  were  developed  to  de terrain e: 

1)  Total  Feature  Worth  (TFW)  to  user 

2)  Percent  of  feature  area  patched 
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3)  Percent  of  feature  with  shattered  slabs 

4)  Equivalent  single  wheel  loads 

5)  Interior  stress  at  the  bottom  of  the  slab 

6)  Cost  of  preparing  pavement  for  overlays 

4. 3. 4.1  Total  Feature  Worth  to  User 

Total  feature  worth  to  the  user  establishes  the  relative  value  of 
the  pavement  feature  to  the  user  with  respect  to  the  entire  network.  It 
is  not  part  of  the  dynamic  programming  analysis  but  it  is  a  function  of 
the  data  used  as  input  into  the  dynamic  programming  analysis.  The  TFV 
is  included  in  the  output  of  the  dynamic  programming  analysis  for  use  in 
the  network  analysis.  Althougn  the  function  TFK  is  used  in  following 
chapters,  an  explanation  is  included  here  because  it  is  output  frco  the 
dynamic  programming  computer  program. 


The  TFW  represents  relative  value  to  the  using  aircraft  of  a 
feature  as  compared  to  another  feature,  whether  on  the  eaoe  base  or 
another  base.  It  is  a  function  of  the  following  variables: 

1)  Pavement  type  (runway,  taxiway,  etc.) 

2)  Pavement  need  (primary,  secondary,  etc.) 

3)  Using  aircraft  type  (bomber,  fighter,  cargo,  etc.) 

4)  Using  aircraft  mission  (alert,  training,  operational, 
etc.) 


i 


76 


5)  Humber  of  using  aircraft  by  type  and  mission. 

Together  these  variables  describe  the  relative  worth  of  one 
particular  pavement  feature  versus  another.  Table  4.1  shows  the 
complete  breakdown  of  each  category  in  addition  to  the  nunerical  value 
assigned  to  each  variable  element.  These  assignments  have  been 
estimated  by  the  author  based  on  6  years  of  experience  in  Air  Force 
pavement  management,  but  should  be  validated  prior  to  implementation. 

The  TFW  of  a  single  feature  is  determined  by  the  following 
relation: 

n 

TFW  *  I  N.  *  (PT-+  PN-  +  AT-+  AM-)  (4.23) 

L=i 

where : 

N.  *  number  of  using  aircraft  by  type  and  mission 
FT^  *  paveaent  type  coefficient  for  using  aircraft 

PH^  =  pavement  need  coefficient  for  using  aircraft 

AT^  =  using  aircraft  type  coefficient 
AM.  =  using  aircraft  mission  coefficient 

i  *  counter  of  different  using  aircraft  types 
n  =  number  of  different  aircraft  types  using  feature 
The  information  necessary  for  the  vortb  calculation  is  readily  available 
at  each  Air  Force  base. 


Table  4,1  K  L  eaten  t.  Weights  of  Totol  Feature  Worth  Calculati 
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*  Relative  weight  of  Characteristic?  as  a  group. 

*#  Relative  weights  of  Criteria  with  in  the  Characteristics. 
Htt  Product  of  Characteristic  and  Criteria  weight. 
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4. 3. 4. 2  Percent  of  Feature  Area  Patched 

Because  the  percent  of  the  feature  patched  is  used  in  the  PCI 
predictive  equation  for  rigid  pavements,  it  was  necessary  to  forecast 
its  value.  Since  it  is  not  a  primary  focus  of  this  research  to  develop 
a  comprehensive  set  of  regression  equations  for  predicting  the  percent 
area  of  feature  patch  as  a  function  of  PCI,  only  a  simple  analysis  using 
91  data  points  from  a  single  airfield  was  used  in  this  research.  These 
points  were  used  to  deteraine  the  coefficients  of  a  third  order 
polynomial : 

PATCH  =  a  *  b*PCI  +  c*PCI2+  d*PCI3  (4.24) 

where: 

PATCH  =  percent  of  slabs  patched  (patch  larger  than  5  ft2) 

PCI  =  pavement  condition  index 

with  the  coefficients: 

a  =  32.204 

b  =  0.1306 

c  =  -0.0098 

d  =  5. 28*10“ 3 

Figure  4.5  depicts  a  graph  of  this  equation. 


79 


Figure  4.5  Percent  o£  Feature 
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4. 3.4. 3  Percent  of  Feature  with  Shattered  Slabs 

The  percent  slabs  replaced  (SR)  in  a  feature  is  also  part  of  the 
PCI  predictive  equation.  For  use  in  the  dynamic  programming  algorithm, 
a  simple  third  order  equation  was  fit  from  data  extrapolated  from  the 
PCI  manual  (27)  as  follows: 

SR  =  a  +  b*PCI  +  c*PCI2+  d*PCI3  (4.25) 

where: 

a  =  88.622 

b  =  -1.993 

c  =  0.012 

d  =  -1 .127*10" 3 

This  equation  predicted  the  number  of  high  severity  shattered  slabs  and 
was  used  as  the  basis  for  determining  the  percent  slabs  replaced  in  the 
feature. 

£.3.4.4  Equivalent  Single  Wheel  Load  (ESWL) 

The  terns  and  ac  in  Equation  3.4  require  the  equivalent  single 
wheel  load  (ESWL)  of  the  using  aircraft.  In  current  practice  the 
procedures  for  determining  ESW’L  use  charts  and  tables(8).  In  order  to 
simplify  the  conputer  program  and  to  save  significant  computation  cine. 


OEPTM  BELOW  PAVEMENT  SUKfACe  fWCHEJ) 


SI 


ESWL.  PtRCtKT  or  LOAD  ON  CONTROLLING  MJMBtR  OF  WMtELS 

O  20  40  60  80  100 


gure  4-6  ESWL  as  Function  of  Aircraft  Type  and  Depth  Below 
Pavenent  Surface  (31) 


£SWU.  PERCENT  Qf  V.OAD  ON  CONTI^tLWC  MOUSES  OF  WHEELS 


0  20  40  «0  80  IOO 


a  third  order  polynomial  was  fit  to  Figures  A. 6  and  A. 7.  The  equation 


is: 

ESHL  =  Cnl  +  Ce2-D  +  Cn3*D2  *  Ca4*D3  (4.^) 

where: 

ESWL  =  equivalent  single  wheel  load  as  a  percent  of  total 
controling  wheel  load  (8). 

D  *  depth  to  point  of  interest  in  pavement  system  in  inches. 
San  =  coefficients  of  polynomial  from  Table  A. 2 
n  =  using  aircraft  index 
n  =  coefficient  index 


A.3.A.:>  Interior  Stress  at  Bottom  of  Slab 


The  term  FAT  in  the  PCI  predictive  Equation  3.2,  is  a  function  of 
the  induced  and  allowable  stresses  ir  the  slab.  It  also  is  current 
practice  to  use  charts  and  graphs  for  determining  the  interior  stress  at 
the  bottom  of  the  slab  (8).  For  ease  in  computation,  instead  of  using 
Figures  A. 8,  A.y  and  A. 10  for  interior  stress  determination,  regression 
equations  were  developed. 


The 

following  equation  form  was  used 

with 

the 

coefficients 

regressed 

for  the  three  cases,  single. 

dual , 

and 

tandum  wheel 

configurations : 
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Table  4.2  Coefficients  for  Equation  4.25  (Cian),  Equivalent 
Single  Wheel  Load 


■  BBVllkl  ■■■■■■  I 

■b*iihibbbbbb| 

■  ■Hill  !  ■■■■■■! 

■bum bib ■■■■■! 


COACrdt  Wenor  S»W 
ond  Slot)  Thctotn  rt 
Grass  Atrcrxtft  ¥tufit 
Single  Geor  Loaned 
at  S*ot>  treeiio* 


Hi 

■  ■■ItUUlHRBBl 

■  sH.ivi'rsaHaal 

ibiiuiiibiiiBBHBHHHH 
■■■lutav  «■■■■■■■■■■■■■■■ 
un  ii  lki  iBimiBiiiBiBBBriri'irir  a  vii 

■  BIBl  lUUlBBBBBBBBBBBBBBPB'JFjB'.Jl'.4BriBr 

■  ■■■  lift  ft'BIBHi|ftfllBBBBBBririV  JfiB^Br.BVJ 
■iifttUftiiaiiQiiBiiiiBrvi^VirflViBii 
■■■■IlitV-VBBBBBBBBBBBBB'jrjrB'JBJB'JBF.BB 
■■iiiiuiiiBBiiBiiaQiir.riirjr>Br.ftrdBri 

■  aaBBlUrilBBBiBBiBBBBB'JfiPBrflrBB^Br.B 

■  BIBB  11  Bft  ■  ■■■■■■■■■riftrfl'iB iBJBMVill 

■  ■■■BlMl  lBBBBBBBBBBBB'iFBriBrjBr.BB’JBB 

■BIBB1 1!  Ill  IBIIBBBMr^irJifflrfiriBMB^ 

■  ■■■BBlkUUlBBBBBBf"r*»jrjBJBJB'JBrflBrJB 


M&wAmwA'ATtTmAWfXMUAnr.mmrmuml 
MHHHiuTirjr«iri>^iBinrur«HBl 
■■■BBiBBiuftiviiiriiri^L/iiMFBjniBaiil 
■■■■■■■■■k^iiLi»rim.ir.a^<Br<Bp^rTvJBHBBBB 
■■■■■■IBBllHU #  WATaWATMT  Il'il^BBBBflg 
tmnmnmunum  H  i  i  tr.wmfAwmwAmum^jmumummM 

MHMB^Va'M'irfl^B'iBfiBB'VBBBBBBDBP 

a’ ArJ  )PJfiBriBMBrilBIBaiHBB| 
BBBBBBBBBBM  «.*4T,k»  l'ririf>irBBr<IBIBBBBBBl| 

■  ■■■■■■■■■/•  ftV'j  riiFJPJBrBBBBBBIIHBl 

■  ■■■■■■■Bfi’i'i'  Til  ,V.W'A9  JBr.flBBBBBBBBBBBP 
IBBBBBpBrf.rr  BA  .  rBTiB^BrBBBBBBBBBBBBHL 

##«•«  u-  .1  >•«!?&■  IT  fl ■  B BBBBBBBBBBB 1 

■  iV  ftr^irBBBBBBBBBBBBBBBl 
BBRBiBB/rr.r  r  a'avssc*  bibbbbbbbbbbbbbbP 


_  i.niL 

BiBirf  r.rjrjr4ra'jm  +  k’^ftii 


BBBBaBBBBBBBBb 


BBBB'iy/f  f  0iMril  <^UVT/v<9BBBBBBBBBBBBB| 

■  Hirfirii^fir.B^B»ii>J^y.^3iBBBBBBBBBBBl 

■  ■■>/, rK^^BBBBBBBBBIBBL 

■  ar/J/f  f  FflinrvBIBrk.:\Vftl^BBBBBIBBIBB| 

■  B/^iiiiirif  i^BIBBBBilkrOmiBBBBI 

ir/y/,ri'iririPir.BiBBBiBftii:<fti^BBili 


■■  iiimin  w  -j.r.rir<r-;,wm'Armr.mr 
lanai'i  *\m\sUuif-  A'iT.a'iv»*A 
lift  liBlIlk/i'i'i 9  A  T*M 'AM '4 Tt  __ 

■■■'■  ■  laiirt/rr  jaA'A9mrjrr,<r-9  jrm'jr.mrmm 
■  am'iii  ■>  M>ii,l'Fir.m'iwfk'£*iTMTAWMrmm-A 
<■■■■■■!■■  inri  I kit.T.t m *r.m-Ac-m VAT *rM-AW , 

llllllll  m  M/ii'i  Ji’rjMr/.tArjNr/r.i’ii 

■  ■■■■■■■lb  '’.’H n  tr.ra'AriyA.rw Ar/^A^m 

li!!!!I!!*  ' '  - 


mm 

■Amumr. 
'jura 

_ „  -BK,mrmm 

SSSSSS!17  r‘  i*’**r‘**r‘r‘rrajw vsraraa'ai* 

. .  vi  >  ri,'«  irritj'irj'iir.KjriHi 

wnmmmww//i  >  <i  s  (ir.i  iii'ir.riM'.iiiaviria 

■  ',i  r.unrvi  I'iMijr.ir.irii'iimi 

'»r& a  At  r,r.rr<iri'1irii'Jjiriau 

■  ■■■■'//rvi  U.',i,'iV»i'<,iri!r  iru  a  a  ana 

■  ■■■'//,  44#  '•J’.ir<'fi'i<-<r.k’ir«r.ir^'{|a»all 

f  j’iraral/iiBiain 

umm’/'j /.r.r.fd  ',uvi',r„'i  jvr4r.ar, •*!■■•■■■■■■ 

*mw/rt  i  trr  lar/ii.fa  4  <a  iriiraraniaaiiiia 

WM'/tJ/SA  JA**  ‘ m  '  r*'M  k  wJBRHaaaHHHHa 

■  r  •  / '/  /  rsi'.r.t  /  r^iTA'A^.a  At*  a*  a  a  asaisiaDiiaaa 

w/vntf.rt  A  /trr  r.rh  iw  a  'arak'Bmiaiiin 
r  r r t > t nr,t t r r r MAVA'  ar.auk/aaiaBBiaia 
/  ,/7/,i'i'ir'/ii'4i,4'ir4.’.0'’iik'v%^iiaaaaB*" 
/  7,  '4  /  r ,  r,  a  /  -a  '4  r4  >"4  b,v:0SS^>  5  a  ■  aa  ■  ■  ■ 

/  /'itta/tfi’i'^r.wA^a^aar^.  ^Vii.'iiaoiaai 
/,  /  'i  '4  4  v  /  r  ft  r.  r.  r  a  a  vb  ta  m  m  a  : ;  ,<1S  w  a  ■  w  •  a  ■  ■ 

//4i’/'i////4'i^a'4anBiiiiaikikii^ak;eiia 

/.r.rI/,>//4'4’4rr4HiuBigig»?Kt7it;(kiii 

»47//-4/»4,ir,a,4aaaaaBBBBaaaaaki>ikat.<»."a 


Figure  4.9 


Interior  Stress  at  the  hottest  of  the  Slab  frost  Dual 
Wheel  Loading  (8) 
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SLBSTRS  -  Cb1*fi  +  Co2*T2  +  C^m  +  C^LGK 

+  Ce5*LGT  +  C^LGL  *  Cn7*IL  +  Cog  (4.27) 


with: 

FI  =  L/T2*Log10(24.165*T°*75/K0‘25) 

F2  =  L/T2 
IT2  «  1/T2 

-  1°810(K) 

LGT  =  log10(T) 

I£L  =  log10(L) 

11  -  1/1 

Cmn  *  regression  coefficients  froo  Table  4.3 
n  =  coefficient  index,  I  thru  8 
n  »  aircraft  wheel  type,  1,  2,  or  3 

where: 

L  =  aircraft  load  in  pounds 

K  *  nod ulus  of  soil  reaction  in  pounds  per  square  inch  per  inch 
T  «  slab  thickness  in  inches 

4.3 .4.6  Cost  of  Preparing  Paveaent  for  Overlay 


The  cost  of  overlaying  a  paveaent  in  the  future  must  reflect  the 
condition  of  the  pavement  feature  at  the  time -of  the  overlay.  The  worse 
condition  a  pavement  is  in  (lower  PCI)  the  more  it  costs  to  repair  the 
paveaent  in  preparation  for  the  overlay.  A  data  set  of  200  features 


Coefficients  for  Equation  4.26  (Cum),  Interior  Stress 
at  the  Bottom  of  tho  Slab 


were  analyzed  from  airfields  around  the  country,  to  determine  the  cost 
of  preparing  the  pavement  feature  for  overlay.  The  following  equation 
resulted  from  regression  of  the  data: 

PREPCST  =  264.70*PCI-1-539  -  0.2506  (4.28) 

where: 

PKEPCST  =  cost  to  prepare  pavement  for  overlay  (1981  $/sy) 

PCI  =  current  pavement  condition  index  of  feature 
between  limits  of  30  and  100 

4.3.3  Dynamic  Programming  Computer  Outputs 

Execution  of  the  computerized  dynamic  programming  model  optimizes 
project  level  major  rehabilitation  activity  plans.  The  model  analyzes 
five  plans  for  each  feature.  The  first  plan  has  no  funds  available  for 
rehabilitation  and  only  routine  maintenance  is  applied.  The  remaining 
four  plans  are  a  functicr  funding  level  limitation.  The  four  funding 
levels  are  equal  intervals  between  whet  it  cobZJ  to  apply  the  cheapest 
activity  (overlay  with  2  inches  of  asphaltic  at  the  beginning 

of  the  analysis  period  and  the  most  expensive  activity  (reconstruction 
with  portland  cement  concrete)  inflated  to  the  end  of  the  analysis 
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The  program  finds  the  particular  set  of  activities,  and  their 
schedules,  that  maximizes  the  area  under  the  utility  weighted  PCI  versus 
time  curve  subject  to  the  major  constraint  of  a  fixed  funding  level. 
Figure  4.11  is  an  example  of  the  output  from  this  program.  Since  the 
output  vas  designed  as  merely  input  into  the  network  analysis  portion  of 
this  research,  it  incorporates  minimal  user  labeling.  Furthermore,  the 
output  contains  only  those  plans  for  each  feature  that  are  unique; 
duplicate  plans  are  eliminated. 

The  output  includes: 

1)  Base  and  feature  name 

2)  Total  benefit  or  performance  of  entire  plan  (area  under 
curve) 

3)  Total  cost  of  plan  in  dollars 

4)  TFW,  total  feature  vorth,  relative  value  of  feature  to  user 

5)  Decision  (selected  activity),  cost  of  decision,  accumulated 
benefit  (area),  and  PCI,  all  as  a  function  of  age 


For  a  typical  airfield,  an  average  of  3  plans  per  feature  are 
generated  for  approximately  150  features.  These  are  optimized  plans, 
that  is,  the  area  under  the  utility  weighted  PCI  versus  tine  plot  has 
been  maximized  with  constraints,  the  most  important  of  which  are 
specified  funding  limitations  over  the  span  of  the  analysis  period. 
Mote  that  all  costs  are  presented  in  present  worth  dollars.  The 
decisions  are  indexed  as  1,  2,  3,  4,  and  5  and  the  same  as  explained  in 
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♦EXAMPLE*  E01A 

BASE  miUIE  BDfEFIT  COST  TFW 


1127. 

0 

48170 

ACE 

DECIS10K 

COST 

8EKEFIT 

PCI 

0. 

1. 

0. 

0 

65 

2. 

1. 

0. 

164 

62 

4. 

1. 

0. 

317 

58 

6. 

1. 

0. 

460. 

54 

8. 

1. 

0. 

592. 

50 

10. 

1. 

0. 

712. 

45 

12. 

1. 

0. 

820. 

41 

14. 

1. 

0. 

916 

26 

16. 

1. 

0. 

9*9. 

31 

18. 

1. 

0. 

1069. 

26 

20. 

1. 

0. 

1127 

21 
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Paragraph  4,3*1. 

On  the  average,  the  computer  program  takes  3  central  processing 
seconds  per  feature  to  generate  the  five  plans  on  the  University  of 
Illinois  Cyber  175  computer. 


4.4  S  unmar 


This  chapter  explained  the  dynamic  programming  procedures  as  they 
applied  to  project  level  optimization  and  discussed  the  mechanics  of  the 
procedures  used  in  this  research.  The  next  chapter  introduces  the 
decision  process  of  long  range  network  level  optimization. 
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CHAPTER  5 

LONG  RANGE  NETWORK  LEVEL  OPTIMIZATION 


5.1  General 


This  chapter  discusses  and  models  the  decision  process  used  by 
pavement  engineers  to  develop  long  range  major  rehabilitative  plans  for 
entire  systems  of  pavement  sections.  The  chapter  includes:  an 
explanation  of  airfield  systans  (networks),  decision  process 
description,  adaption  of  zero-one  integer  linear  programming  to  the 
decision  process,  a  specification  of  the  optimization  criterion  and 
optimization  constraints,  and  an  explanation  of  the  specific  integer 
programming  method  used  to  obtain  the  optimal  decisions. 
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5*2  Airfield  Pavement  Systems  (Hctworks) 


HetworV  level  planning  involves  considering  several  long  range 
rehabilitation  projects  for  each  feature  on  the  airfield.  Although  eac* 
of  these  projects  may  be  optimal  for  a  particular  feature,  there  stil? 
remains  the  problem  of  selecting  the  best  of  these  projects  to  satisfy 
the  overall  needs  of  the  using  aircraft  as  touch  as  possible  vithin  the 
bounds  ot  a  limited  budget  for  the  entire  network.  Using  long  range 
network  level  optimization,  airfield  pavement  engineers  can  select  the 
best  set  ot  long-range  projects  for  the  network. 


5.2.1  Airfield  Pavement  Ketvork  Definition 


Definition  of  an  airfield  pavement  network  is  common  to  ooth 
military  and  civilian  airfields.  An  airfield  pavement  network  is  a 
grouping  of  two  or  more  airfield  pavement  features.  A  network  can 
consist  of  only  a  portion  of  an  airfield,  i.e.,  runway  or  an  spron,  or 
an  entire  airfield,  complete  with  all  the  features  including  all  the 
aprons,  taxivays,  and  runways.  Still  a  third  definition  group!  several 
distinct  airfields  many  miles  apart  into  a  network*  These  network 
groups  can  be  located  statewide,  command  vide,  country  vide  or  scattered 
throughout  the  world. 
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For  tbe  most  part,  a  network  level  analysis  encompasses  only  t**ose 
pavement  features  under  the  control  of  a  specific  decision  authority. 
One  airfield  would  generally  not  conduct  a  network  analysis  for  itself 
combined  with  another  airf ie* i  not  under  their  control.  Yet  on  the 
other  hand,  a  governing  agency  responsible  for  several  airfields  might 
justifiably  include  all  the  features  on  all  the  airfields  within  their 
jurisdiction  into  their  network  analysis. 


5.2.2  Air  Force  Network  Definitions 


The  U.S.  Air  Force  command  structure  is  logically  set  up  for 
networks  defined  in  several  ways.  First,  individual  airfields 
constitute  the  most  prevalent  network.  Next,  the  Air  Force  is  divided 
in  separate  operating  commands  with  jurisdiction  over  airfields  spread 
throughout  the  world.  These  command  divisions  could  delineate  a  second 
level  ot  networks  in  the  Air  Force.  Lastly,  Air  Force  headquarters 
reserves  the  highest  level  of  funding  approval  for  itself  as  it  directly 
petitions  for  funds  from  the  Congress.  Therefore,  anotner  definition  of 
a  network  could  be  all  airfields  in  the  Air  Force. 

Summarizing,  the  Air  Force  has  three  basic  network  levels  to  which 
optimization  techniques  could  be  applied: 

1)  base 

2)  command 

3)  Air  Force  wide 
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5-3  Network  Level  Decision  Procees 

Airfield  pavement  engineers  are  faced  with  coupler  decision  process 
when  developing  long  range  network  level  major  rehabilitation  plans-  At 
this  level,  they  must  select  from  a  list,  those  projects  that  best 
satisfy  the  needs  of  the  user  within  the  bounds  (constraints)  placed 
upon  them  (usually  by  the  user).  These  might  include,  but  are  not 
limited  to: 

1)  budget  limitations  (included) 

2)  relative  value  of  features  to  user  (included) 

3)  quality  of  rehabilitative  activity  to  feature  (included) 

4)  political  influence  on  affected  features 

5)  operational  impact  on  timing  of  rehabilitation  activity 

5-3-1  Decision  Process  Description 

At  the  network  level,  the  best  projects  must  be  selected  within  a 
given  funding  limit-  One  method  for  selecting  the  best  projects  can  he 
determined  by  prioritizing  the  features  according  to  user  needs-  With 
only  one  project  per  feature  (determined  by  a  decision  tree)  the 
projects  can  be  ranked  according  to  user  needs  using  the  feature 
prioritizing-  Budgets  can  then  be  allocated  down  the  list  until  funds 
are  depleted-  However,  this  decision  process  does  not  consider: 


1 
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1)  Relative  usefulness  to  user  between  features  in  network. 

2)  Relative  structural  benefit  of  activities  to  a  feature 
between  projects  of  the  same  feature  and  between  projects  of 
different  features. 


Pavement  engineers  have  difficulty  in  manually  selecting  the  best 
set  of  projects  given  a  limited  amount  of  money,  because  there  are  too 
many  possible  combinations  of  projects  for  them  to  analyze.  Each 
combination  provides  a  corresponding  level  of  satisfaction  to  the  user 
and  to  the  pavement  structure.  To  best  spend  the  limited  funds,  the 
selected  projects  should  maximize  (optimize)  the  level  of  satisfaction 
for  both  of  these  needs  (user  and  structure).  Using  OR  techniques,  the 
set  of  projects  can  be  selected  so  that  the  highest  possible 
satisfaction  is  achieved. 


5.3.2  Hetvork  Level  Decision  Variables 


There  is  only  one  type  of  decision  variable  in  the  model  used  to 
solve  the  network  level  optimization  problem:  there  is  one  variable  of 
this  type  for  each  alternative  project  for  each  feature.  In  the  case  of 
150  features  averaging  three  alternative  projects  per  feature,  there 
would  be  450  decision  variables  in  the  problem.  These  variables  are 
constrained  by  definition  to  take  on  the  value  of  either  zero  or  one. 
If  a  decision  variable  takes  on  the  value  one,  then  the  project  is 
selected.  If  it  is  zero,  then  it  is  not  selected. 
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5.3.3  Hetwork  Level  Optimization  Criterion 


The  network  level  procedures  developed  in  this  research  reflect 
long  range  planning  because  they  select  among  projects  already  optimized 
over  a  long  analysis  period.  The  network  level  decision  process  does 
not  differentiate  time  as  a  variable  in  the  analysis.  If  the  feasible 
projects  reflected  activities  scheduled  for  the  current  year,  then  the 
network  optimization  would  also  only  determine  what  to  do  for  the 
current  year.  The  network  analysis  procedure  meets  the  future  needs 
only  if  the  projects  subject  to  selection,  are  long  range  plans. 

The  network  optimization  ciiterion  is  a  function  of  two  parameters: 

1)  Importance  of  feature  to  the  user 

2)  Benefit  of  the  project  to  the  structural  performance  of 
the  feature. 

5.3 .3.1  Importance  of  Feature 

The  total  feature  worth  (TFW)  is  the  importance  or  value  placed 
upon  a  feature  by  the  user.  TFW  is  a  function  of:  pavement  type, 
pavement  need,  using  aircraft  type,  using  aircraft  mission,  and  the 
number  of  aircraft  by  type  and  mission  (previously  defined  in  Chapter 

3).  This  parameter  quantifies  the  value  of  one  feature  relative  to 
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another.  Features  need  not  be  adjacent  or  even  located  on  the  same 
airfield.  TFW  is  more  than  just  a  prioritized  ranking  of  the  most 
important  feature  to  the  least  important  feature.  Its  magnitude,  as 
compared  to  each  other,  signifies  the  relative  value  of  each  feature  to 
the  user.  If  the  TFW  of  a  feature  is  tvice  that  of  another,  then  the 
value  of  the  first  feature  to  the  user  is  tvice  that  of  the  other. 

TFW  is  computed  for  each  feature  considering  the  same  data  used  for 
the  development  of  long  range  rehabilitative  plans  at  the  project  level* 
The  actual  TFW  value  for  each  feature  is  the  subjective  opinion  of  the 
author  based  on  bis  six  years  experience  as  an  Air  Force  pavement 
engineer. 


5*3. 3. 2  Structural  Benefit  to  Feature 


The  structural  benefit  of  the  project  to  the  feature  (performance) 
is  defined  as  the  area  under  the  utility  weighted  PCI  versus  time  curve. 
This  parameter,  generated  from  the  project  level  optimization  process, 
represents  the  structural  benefit  to  a  feature  from  a  specific  set  of 
rehabilitated  activities*  It  is  not  merely  a  prioritized  ranking  of 
projects  from  the  same  feature,  but  it  indicates  the  relative 
performance  of  the  scheduled  activities  of  projects  fra*  not  only  the 
same  feature  but  also  different  features. 


5. 3. 3. 3  Optimization  Criterion 

The  optimization  criterion  in  this  research  uses  a  combination  of 
TFW  and  structural  performance.  Using  either  alone  vould  neglect 
important  considerations  measured  by  the  other  parameter.  Network  plans 
optimized  with  only  TFW  vould  not  consider  structural  consequences, 
either  positive  or  negative,  with  different  levels  of  funding. 
Additionally,  long  range  network  plans  optimized  vith  just  structural 
performance  vould  not  consider  the  needs  of  the  user  as  it  selected 
projects.  Conceivably,  if  structural  performance  alone  were  used, 
rarely  used  pavement  features  could  receive  more  activity  than  an 
important  highly  used  feature. 

The  optimization  criterion,  used  in  this  research,  is  the  TFW 
veighted  structural  performance.  For  each  project  in  the  network 
analysis,  the  product  of  TFW  and  structural  benefit  is  computed  and 
represents  its  contribution,  if  selected,  to  the  overall  objective. 

5.3.4  Network  Level  Optimization  Constraints 

For  this  research,  tvo  types  of  constraints  vere  included  in  the 
network  level  optimization  model.  First  and  most  isportantly,  the 
selected  plans  must  not  exceed  a  pre-established  funding  limit.  And 
secondly,  only  one  project  per  feature  should  be  included  in  each 
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selected  plan.  Operational  constraints  (i.e.,  pre-specif ied  features 
must  either  be  selected  or  denied  activities  as  a  unit)  could  easily  be 
added,  but  are  not  included  here. 


5.4  Zero-One  Inteeer  Linear  Pro 


The  decision  variables  in  ©any  problems  must  sometimes  be  limited 
to  integer  values.  For  example,  people,  machines,  animals,  cars, 
trucks,  etc.  can  only  be  allocated  in  vhole  (integer)  quantities. 
Integer  optimization  problems  vith  only  linear  functions  in  the  model 
are  classified  as  integer  linear  programming  problems.  Further 
limitations  can  arise  if  the  decision  variables  r epresent  yes  or  no 
decisions.  For  example,  such  a  decision  variable  could  represent  the 
selection  of  a  specific  policy.  The  decision  variable  is  not  only 
restricted  to  being  an  integer,  but  it  is  also  restricted  to  only  two 
values,  zero  or  one  (zero  means  no  and  one  means  yes).  This  type  of 
integer  programming  problems  vith  zero-cne  decision  variables  and  linear 
functions  is  called  zero-one  integer  linear  programming  (20). 

In  this  research,  the  network  level  decision  process  is  modeled  as 
a  zero-one  integer  linear  programming  problem.  The  decision  variables 
represent  the  projects  for  all  the  features  in  the  network.  The  values 
for  each  decision  variable  can  only  be  zero  or  one,  where  zero  means  the 
project  was  not  selected  and  one  means  the  project  was  selected. 
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Modeling  the  network  level  decision  process  in  this  manner  permits 
the  selection  of  projects  that  maximizes  the  objective  function 
(performance  criteria)  while  maintaining  a  pre-established  budget. 
However,  obtaining  an  optimal  solution  by  integer  programming  techniques 
is  infeasible  because  most  problems  will  contain  too  many  decision 
variables  to  be  solved  in  a  reasonable  amount  of  computer  tine, 
therefore,  a  heuristic  solution  process  has  been  used  to  obtain 
solutions  which  are  generally  close  to  optimal. 


5.5  Zero-One  Integer 


.  of.  Metwork  Decision 


IlSKga 


The  zero-one  integer  linear  programming  model  (here  after  referred 
to  as  the  IP  model)  is  composed  of  the  objective  function  and 
constraints.  These  two  parts  are  mathematical  expressions  of  the 
decision  variables  (long  long  range  major  rehabilitative  projects  to 
undertake). 


5.5.1  Objective  Function 


The  objective  function  of  the  network  level  decision  process  can  be 
modeled  as: 

Maximize: 


Z  Z 
i  i 


* 


*  TFWj 


(5.1) 


S-:  .1  £'— 


decision  variable  of  the  jth  project  and  ith  feature 
performance  of  the  jth  project  and  ith  feature 
(area  under  utility  veighted  PCI  curve) 
total  feature  worth  of  the  ith  feature 
feature  counter 


project  counter 


The  product  of  the  parameters  and  TF#£  represent  both  the 
structural  performance  of  the  activities  in  the  project  to  the  feature 
and  the  relative  need  of  the  project's  feature  to  the  user.  P^. 

TFW^  are  both  determined  during  the  project  level  analysis  as  explained 
in  Chapter  3. 


5.5.4  Constraint  Functions 


Three  constraints  are  included  in  the  model  used  in  this  research. 
First  the  total  cost  of  the  selected  projects  must  cot  exceed  a 
pre-established  limit.  This  can  be  mathematically  expressed  as  follows: 


l  IX..  *  C;;  <  MAXCOST 


decision  variable  of  the  jth  project  and  ith  feature 
cost  of  the  jth  project  and  ith  feature 
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For  all  projects  selected  (X—  **  1)>  this  sunmation  vill  determine  the 

total  network  cost. 


Secondly,  it  is  necessary  to  mathematically  restrict  selection  to 
no  more  than  one  project  per  feature.  This  can  be  accomplished  with: 


£  £  i  for  Aii  i  (5.3) 

j 


where: 

Xij  “  decision  variable  of  the  jth  project  and  ith  feature 
This  expression  limits  the  selection  of  no  more  than  one  project  per 
feature,  but  the  functions'  inequality  does  not  force  the  model  to 
select  a  project  for  every  feature. 


Lastly,  the  requirement  that  the  decision  variables  are  limited  to 
the  integer  values  of  zero  or  one,  mathematically  restricts  the  optimal 
solution  to  whole  projects. 


5.6  I?  Model  Solution 


Because  of  the  limitation  of  the  computational  capabilities  of 
current  integer  programming  methods,  the  zero-one  integer  model  in  this 
research  was  solved  using  a  highly  efficient  heuristic  (36).  Resource 
allocation  problems  are  commonly  modeled  with  zero-one  integer 
programming.  Solving  zero-one  integer  models  is  difficult  with  a  large 


106 


number  of  decision  variables  because  such  problems  can  not  be  solved  in 
a  reasonable  amount  of  time  (34,20,36,37)*  Since  the  network  level 
allocation  problems  developed  in  this  research  have  a  large  number  of 
decision  variables,  it  is  infeasible  to  solve  them  optimally. 
Heuristics  have  been  developed  which  do  not  guarantee  an  optimal 
solution,  but  which  provide  a  solution  reasonably  close  to  optimal. 


5.6.1  Toyoaa  Algorithm 


Toyoda"s  heuristic  for  obtaining  approximate  solutions  to  zero-one 
programming  problems  is  highly  successful  in  determining  optimal 
solutions  to  specific  classes  of  problems  (38,39).  In  general,  zero-one 
integer  programming  problems  solvable  with  this  algorithm  can  be 
described  as  follows: 


a 

MAXIMIZE:  Z  -  Z  R.  *  x-  (5.4) 

L=1  1  1 

m 

SUBJECT  TO:  v  H..  *  %i  <  Li  for  j  *  l,2,...,n  (5.5) 

L=1 

AND:  X.  «  0  or  1  for  i  =  l,2,...,m  (5.6) 


where: 

Z  **  objective  function  value 
*  performance  of  project  i 
X.  =  decision  variable  r presenting  project  i 


n 


i 

i 

with: 


magnitude  of  resource  j  required  for  project  i 

upper  limit  of  resource  j 

number  of  restricted  resources 

number  of  projects 

index  of  projects 

index  of  resources 


when: 


>  0 


Hij  >° 
Lj  >  0 


Xj  «  1  then  project  is  selected 
X-  •  o  then  project  is  not  selected 


A  later  section  contains  an  explanation  of  the  mechanics  of  the 
algorithm  and  its  performance  in  determining  the  optimal  solution. 


The  network  level  project  selection  model  is  within  the  scope  of 
the  Toyoaa  algorithm.  Using  it,  a  set  of  projects  can  be  identified 
that  maximizes  the  satisfaction  of  the  user  and  the  needs  of  the 


pavement  structure. 
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5.6.2  Computer  Program  Description 


A  computer  program  vas  written  for  using  Toyoda's  algorithm  in 
selecting  the  optimal  t>et  of  projects  for  a  network  that  maximizes  the 
product  of  the  project's  structural  performance  (area)  and  the  effected 
feature'6  value  to  the  user  (TFW).  The  program  was  written  in  CDC's 
Fortran  V  with  826  lines  and  15  subroutines  supporting  the  main  program. 
Figure  5.1  is  a  program  flow  chart  showing  the  major  elements  of  the 
Toyoda  algorithm  for  solving  large  IP  problems  of  the  appropriate  form. 

After  the  problem  data  are  entered  into  the  program  a  candidate 
project  list  is  established.  If  candidate  projects  are  present  then 
effective  gradients  (defined  in  next  section)  for  all  the  projects  are 
computed.  The  project  with  the  maximum  effective  gradient  is  selected 
as  part  of  the  optimal  solution.  Model  constraints  are  checked,  and  a 
new  candidate  project  list  is  established.  If  there  are  projects  left 
to  consider  then  the  gradients  are  again  determined  for  each  project, 
the  project  with  maximum  gradient  is  selected,  etc.  until  there  are  no 
candidate  projects  remaining.  The  results  of  the  algorithm  are  the 
selected  projects  from  each  interaction. 


i 


no 


5.6*3  Algorithm  Explanation 


Toyoda's  algorithm  is  a  simple  and  quick  method  for  obtaining 
approximate  solutions  to  large  scale  zero-one  programming  problems  of  a 
certain  form.  The  method  does  not  use  enumeration,  as  do  most  other 
methods,  but  instead  uses  a  measure  of  preferability  (effective 
gradient)  to  change  the  decision  variables  from  zero  to  one  (selected 
projects)  (39). 

The  preferability  of  each  feasible  project  is  the  effective 
gradient  or  the  unit  gain  (as  measured  by  the  objective  function)  per 
the  resources  expended.  When  the  problem  contains  only  one  resource  it 
is  easy  to  find  the  preferability  of  each  project.  When  more  than  one 
resource  (model  constraint)  exists,  it  is  not  easy  to  evaluate  and 
compare  between  projects  the  unit  gain  in  the  objective  per  resources 
expended.  Toyoda's  algorithm  uses  a  penalty  vector  to  establish  a 
scalar  quantity  representing  the  gain  in  objective  per  all  the  resources 
expended.  This  effective  gradient  is  then  used  to  select  a  project  for 
the  final  solution.  The  effective  gradient  for  the  remaining  projects 
changes  as  projects  are  selected  one  at  a  time. 
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5.6. A  Computer  Program  Inputs 


Data  for  the  computer  program  can  be  input  in  three  vays: 
interactively  by  the  user,  a  previously  interactively  generated  data 
file,  and  a  condensed  data  file  generated  by  the  computer  based  on 
output  from  the  project  level  optimizations. 


5.6 .4.1  Interactive  Dser  Input 


The  program  can  optionally  interactively  praapt  the  user  for  all 
the  information  to  execute  the  algorithm.  For  small  problems  (small 
networks)  this  method  is  satisfactory.  For  large  problems,  interactive 
active  data  input  is  time  consuming  and  tedious. 

The  program  prompts  for  the  following  information: 

1)  Problem  title 

2)  Humber  of  projects 

3)  Humber  of  resources  (constraints) 

4)  Project  requirements  of  each  resource 

5)  Resource  limits 

6)  Project  response  (objective  function  contribution) 

7)  Output  format  instructions 


5.6.4.Z  Interactively  Generated  Data  File  Input 

As  the  user  responds  to  the  prompts  during  an  interactive  input 
session,  the  program  automatically  generates  another  computer  data  file 
with  all  the  input  data.  At  a  later  session,  the  user  can  prompt  the 
computer  to  read  all  the  input  data  froo  this  previously  generated  file. 
In  this  manner  the  user  can  save  time  and  edit  the  data  file  for 
mistakes  and/or  make  changes  for  another  analysis. 

5.6. 4.3  Computer  Generated  Condensed  Data  File  Input 

The  L „o  previous  methods  of  data  input  into  the  Toyoda  algorithm 
are  time  consuming  and  tedious,  for  even  normal  sized  networks.  To 
input  normal  and  larger  sized  network  problems  efficiently,  a  computer 
program  was  written  that  reads  the  data  from  a  file  generated  by  the 
project  level  dynamic  programming  algorithm  and  creates  a  condensed  data 
set  for  the  Toyoda  program.  All  input  parameters  for  the  Toyoda  program 
are  established  except  the  network  funding  limit.  The  network  funding 
limit  is  interactively  input ed  by  the  user  to  preclude  running  the 
condensing  program  repeatedly  for  analysis  of  the  same  network  at  many 
different  funding  levels. 

The  condensing  program  if,  written  in  CDC's  Fortran  V  and  has  384 
lines  with  5  subroutines  supporting  the  main  program.  Figure  5.2  is  a 
flowchart  depicting  the  main  functions  of  the  computer  program. 


INITIALIZE 

PROGRAM 


Figure  5*2  Computer  Program  Flowchart  for  Condensing  Project  Level 
Dynamic  Programming  Data  fo~  Toyoda  Algorithm 
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5*6*5  Toyoda  Computer  Prograa  Outputs 


The  Toyoda  computer  prograa  identifies  those  projects  that  optimize 
the  project's  structural  performance  and  affected  features'  value  for 
the  entire  network  subject  to  a  pre-established  funding  limit  in  present 
worth  dollars*  The  same  condensed  data  can  be  rerun  at  different 
funding  levels  established  by  the  user*  Figure  5.3  is  an  example  of  the 
output  from  the  Toyoda  algorithm  program.  Since  this  information  is 
used  as  input  into  a  report  generating  program  it  currently  has  minimal 
user  labeling.  The  output  contains  the  following  information: 

1)  Problem  title  with  optimization  criteria 

2)  Listing  of  selected  projects 

3)  Objective  function  value  of  selected  projects  only 

4)  Resource  1  limitation  (network  funding  limit) 

5)  Resource  1  used  (amount  of  network  funding  limit  expended 
on  this  set  of  projects) 

The  index  number  identif ing  each  project  is  assigned  in  the  condensing 
program  and  is  interpreted  in  the  report  generating  program  discussed 
next. 


The  condensing  program  requires  little  computer  time  to  execute. 
For  a  network  with  approximately  300  projects  on  100  features,  the 
program  takes  between  3  and  4  central  processing  seconds  to  generate  the 
condensed  data  set  on  the  University  of  Illinois  Cyber  175  computer. 


.5 


EXAMPLE  PROBLEH  WITH  OPTIMIZED  WITH  TFV*PERFOBHANCE 

SELECTED  PROJECTS 
1 
5 
7 

11 

13 

MAXIMIZED  OBJECTIVE  FUNCTION  VALUE 
429306570. 

RESOURCE  1  LIMITATION 
750000. 

RESOURCE  1  USED 
366428. 


Figure  5.3  Output  Exar.ple  for  Toyoda  Algorithn 
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The  Toyoca  algorithm  requires  core  tine  and  it  is  e  function  of  the 
pre-established  funding  limit.  Figure  5.4  depicts  the  central  process 
time  as  a  function  of  funding  limit. 

5.7  Hetvorlc  Analysis  Report 


A  computer  program  was  written  to  generate  a  series  of  three 
reports  for  each  level  of  network  funding  based  on  the  results  of  the 
computer  programs  for  the  project  and  network  level  analyses.  The 

program  reads  from  two  files  (one  each  from  the  dynamic  programming  and 

Toyoda  programs)  and  compiles  the  results  into  a  useable  format.  For 
each  network  funding  limit  the  following  reports  are  compiled: 

1)  Selected  Project  Listing 

2)  Features  Without  Projects 

3)  Network  Summary 

The  report  program  is  written  in  CDC's  Fortran  V  and  consists  of  838 
lines  with  12  subroutines  supporting  the  main  program. 

5.7.1  Selected  Project  Listing 


Figure  5.5  is  an  example  of  the  Selected  Project  Listing  (SPL) 
report  for  the  sane  problem  discussed  in  a  later  section.  For  each 


Central  Proceaaing  Tine  on  University  of  Illinois  Cyber 
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Figure  5.5  Selected  Project  Listing  ExsopL 


level  ol  network  funding  in  the  analysis,  an  SPL  is  compiled.  It  lists 
all  the  projects  that  the  Toyoda  algorithm  selected.  Each  project  is 
listed  with  the  following  data: 


1)  Base  and  feature  name 

2)  Total  cost  of  project  in  present  worth  dollars 

3)  Benefit,  defined  as  the  structural  performance  of  the 
activities  during  the  analysis  period  (total  area  under  the 
utility  weighted  PCI  versus  tine  plot) 

4)  Worth,  defined  as  the  relative  value  of  the  affected  feature  to 
the  user  (function  of  the  aircraft  type,  mission  and  number, 
plus  the  pavement  type  and  need). 

5)  PCI,  pavement  condition  index  as  a  function  of  time 

6)  Decision,  optimized  activities  spanning  the  analysis  period 

7)  Cost,  expenditures  for  activities  during  the  analysis  period  in 
present  worth  dollars 


Preceding  the  project  listings,  the  following  network  analysis  data 
is  anotated: 

1)  Report  title 

2)  Decision  legend 

3)  Hetwork  description  with  optimization  criteria 

4)  Hetwork  funding  limit 

5)  Aoount  spent  in  this  listing  on  present  worth  dollars 

6)  Objective  function,  value  of  the  total  network  plan. 
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5.7.2  Features  Without  Projects  Listing 


Figure  5.6  is  an  example  of  the  Features  Without  Projects  (FWP) 
report  for  the  sample  problem.  The  FWP  report  lists  all  the  features  in 
the  network  that  did  not  receive  any  work  from  the  selected  piojects  and 
the  consequences  of  these  decisions  on  the  specific  features.  Each 
feature  is  listed  with  the  following  data: 

1)  Base  and  feature  name 

2)  Benefit,  defined  as  the  structural  performance  of  the 
feature  with  only  routine  maintenance  during  the  analysis 
period  (total  area  under  the  utility  weighted  PCI  versus  time 
curve) 

3)  Rel  Value,  same  as  worth,  defined  as  the  relative  value  of  the 
feature  to  the  user  (function  of  the  aircraft  type,  mission, 
and  number,  plus  the  pavement  type  and  need) 

4)  PCI,  pavement  condition  index  as  a  function  of  time 

Note  as  with  the  previous  report,  the  listing  is  preceded  by  the  network 

description,  funding  limit  and  amount  spent.  As  anotated  on  the  report, 

those  features  whose  PCI  falls  below  40  are  astericked  and  flagged, 

emphasizing  the  serious  condition  of  these  features  to  the  using 

pavement  engineers.  As  with  the  SPL,  a  FWP  is  also  generated  at  each 
network  funding  level. 
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5.7.3  Network  Summary  Report 

Figure  5.7  is  an  example  of  Network  Summary  (NS)  report.  The  NS 
report  is  a  summary  of  data  for  a  specific  funding  level  of  the  network 
analysis.  Again  it  identifies  the  report  with  a  network  description, 
funding  limit  and  amount  spent.  Then  the  following  information  is 
compiled  as  a  function  of  time  during  the  analysis  period: 

1)  PCI,  an  average  pavement  condition  index  weighted  with 
total  feature  worth  (TFW) 

2)  Cost,  amount  spent  on  the  network  in  present  worth  dollars 

3)  Accumulated  Cost,  the  summed  total  amount  spent  on  the  network 
in  present  worth  dollars 

4)  New  Features  w/PCI  40,  the  number  of  additional  that  have  just 
fallen  below  a  PCI  of  40 

5)  Poor  Features  (FCI  40),  the  total  number  of  features  in  the 
network  with  a  PCI  less  than  40 

6)  Benefit,  total  benefit  from  all  the  features  in  the  network 
during  that  time  interval 

7)  Accumulated  Benefit,  the  summed  total  benefit  from  all  the 
features  in  the  network 


Note,  Benefit  is  defined  as  the  area  under  the  utility  weighted  PCI 
versus  time  plot  for  each  feature  over  the  entire  analysis  period  or 
during  a  specific  interval  of  time. 


Figure  5.7  Network  Summary  Example 


The  report  program  generates  all  three  listings  for  each  level  of 
network  funding.  It  takes  approximately  7  central  processing  seconds 
for  a  single  funding  level  report  on  the  University  of  Illinois  Cyber 
175  computer. 


5.8  Sample  Airfield  , 


As  an  illustration,  a  small  sample  airfield  was  analyzed  with  these 
procedures.  The  sample  airfield  has  only  ten  features  so  the  size  of 
the  problem  do  no.  hinder  understanding  the  decision  process  (Figure 
5.8).  The  inputs  for  each  feature  were  extracted  from  real  Air  Force 
base  features  subjected  to  light  load  aircraft.  All  the  feature  data 
was  run  through  the  dynamic  programming  algorithm  and  these  results  were 
input  into  the  Toyoda  network  analysis  programs.  All  reports  were 
generated  for  funding  levels  between  zero  and  3  million  present  worth 
dollars  (at  §250,000  intervals).  The  resulting  objective  function  (TFW 
weighted  performance)  values  are  plotted  on  Figure  5.9.  As  a 
comparison,  the  same  sample  airfield  was  managed  manually  as  follows. 
Using  three  different  condition  levels  (PCIs  of  40  ,  60,  and  80)  a 
separate  network  analysis  was  completed.  An  activity  (2  inch  asphalt 
overlay)  was  scheluled  for  any  feature  on  the  sample  airfield  when  its 
condition  reached  the  pre-established  minimum  condition  level  (40,60,80 


PCI).  The  analyses  were  carried  out  for  20  year*  (the  same  as  for  the 
Toyoda  analyses).  For  each  of  the  three  analyses,  the  objective 
function  value  was  calculated  and  plotted  on  Figure  5.9. 


.25 


.gure 


Figure  5.9  Objective  Function  Versus  Network  Funds  Spent  for 
Example  Airfield 
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Note,  that  in  each  case,  for  the  sane  aaount  of  money  spent,  the 
methods  developed  in  this  research  nearly  double  the  objective  function 
values  obtained  manually.  Or,  from  another  perspective,  for  any  level 
of  the  objective  function,  the  cost  of  the  optimally  selected  projects 
were  less  than  half  the  cost  of  the  manually  selected  methods  to  achieve 
the  same  objective  function  value. 

Examination  and  comparison  of  the  results  (which  feature  activities 
were  scheduled  at  various  funding  levels)  for  both  methods  revealed 
several  reasons  for  the  vast  differences  illustrated  in  Figure  S.9. 

1)  Activity  assignment  vas  optimized  at  the  project  level  for 
the  method  derived  in  this  research. 

2)  Manual  network  analysis  selected  projects  irrgardless  of  the 
relative  value  of  feature  to  the  user. 

3)  Manual  analysis  also  does  not  take  into  consideration  the 
structural  benefit  (performance)  of  a  selected  activity 
schedule  (project). 


This  small  example  illustrates  the  value  of  the  programs  developed 
in  this  research:  they  provide  substantially  better  ways  of  spending  the 
same  money  and/or  maintaining  an  established  condition  for  a  reduced 
amount  of  money. 
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5.9  Stnmarv 

This  chapter  introduced  the  network  level  decision  process  by 
explaining  airfield  pavement  networks,  the  decision  process,  use  of 
zero-one  integer  programming,  the  optimization/constraint  criteria  for 
the  network  decision  process,  and  solution  by  Toyoda's  heuristic.  A 
general  discussion  of  the  cocaput er  programs  used  in  the  netw>v  '  analysis 
included:  Toyoda's  algorithm,  data  condensing,  and  report  generating. 
The  next  chapter  discusses  an  application  of  the  project  and  network 
level  analysis  procedures  to  McClel lan  Air  Force  Base. 
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CHAPTER  6 

RESEARCH  APPLICATION 


6.1  General 


In  this  chapter,  the  procedures  described  in  the  preceding  chapters 
are  applied  to  McClellan  Air  Force  Base  (HArB),  California.  The  entire 
process  is  presented  froo,  the  initial  data  collection  through  an 
analysis  of  the  final  reports* 


McClellan  AFB  is  located  eight  miles  northeast  of  downtown 
Sacranento,  California  in  the  northern  portion  of  the  Great  Valley  of 
California  (Sacranento  Valley)  and  near  the  western  edge  of  the  Sierra 
Nevada  foothills*  The  airfield  pavements  complex  consists  of  a 


north-south  runway,  parallel  taxiway  and  eight  connecting  taxiways  on 
the  east  side,  five  parking  and  maintenance  aprons  on  the  east  side, 
five  taxiways  and  three  aprons  on  the  west  side,  and  several  maintenance 
hanger/ shelter  access  rasps  located  around  the  airfield*  Figure  6*1  is 
a  layout  of  the  airfield  with  all  the  features  identified  (40)* 

6.2  Data  Collection 

Input  data  for  the  analysis  were  collected  from  several  sources: 

1)  Airfield  pavement  evaluation  report  (40) * 

2)  Pavement  condition  information  from  U.S.  Air  Force 

Logistics  Command  Engineering  and  Services  Office  at 

Wright -Pat ter son  AFB,  Dayton,  Ohio. 

3)  Conversations  with  Air  Force  Logistics  Command  pavement 
engineer  and  McClellan  AFB  pavement  engineer. 

4)  Airfield  layout  data  from  Air  Force  Engineering  and  Services 
Center,  Tyndall  AFB,  Panama  City,  Florida. 

Data  were  collected  via  telephone  conversation,  mail,  and  personal 
visit.  The  pavaent  condition  information  was  obtained  from  a  complete 
PCI  survey  of  all  the  airfield  features,  stored  on  the  command  computer 
data  system  at  Dayton,  Ohio;  and  the  airfield  pavement  evaluation 
information  was  obtained  from  reports  located  at  Dayton,  Ohio;  Panama 
City,  Florida;  and  the  base,  Sacramento,  California. 


Airfield  Layout  Plan  of  McClellan  APB,  California  (40) 
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The  necessary  data  were  extrapolated  free  these  sources  and  coded 
for  the  University  of  Illinois  Cyber  175  computer.  Figure  6.2  is  a 
sample  of  the  coded  data  source  file  for  McClellan  AFB.  The  data  by 
line  are: 

1)  Base  and  feature  name,  pavement  type. 

2)  Feature  length,  width  and  age. 

3)  If  rigid  then  slab  size,  modulus  of  rupture,  modulus  of 
soil  reaction,  slab  thickness,  asphalt  overlay  thickness, 
current  PCI,  percent  slabs  replaced,  freezing  index,  percent 
slabs  patched,  and  average  annual  taperature. 

4)  If  flexible  the,  number  of  layers  in  pavement,  asphalt 
thickness,  age  of  overlay. 

5)  If  flexible  then,  layer  type  code,  layer  thickness,  and  layer 
CBR.  (a  line  for  each  layer) 

6)  If  flexible,  current  PCI,  average  annual  temperature,  freezing 
index. 

7)  Payment  use,  pavement  need. 

8)  Maximum  asphalt  overlay  thickness,  reconstructed  pcc  thickness, 

o 

overlay  cost  per  yd  per  inch  thick,  reconstruction  cost  per 
yd2,  reconstructed  modulus  of  rupture,  reconstructed  modulus  of 
soil  reaction,  reconstructed  slab  size. 

After  the  feature  property  data  have  been  coded  into  the  computer, 
the  dynamic  programming  algorithm  is  used  to  generate  optimized  long 
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Figure  6.2  Abbreviated  Data  Source  File  for  McClellan  AFB 


range  major  rehabilitation  plans  for  each  feature.  The  algorithm  ran 
for  340  central  processing  seconds  on  the  University  of  Illinois  Cyber 
175  computer.  Figure  6.3  is  an  example  of  the  output  from  this  program 
for  McClellan  AFB.  For  each  feature  five  funding  levels  were  analyzed. 
As  shown  for  feature  R01A,  on  the  figure,  the  first  project  for  each 
feature  is  always  limited  to  zero  funding.  Four  other  funding  levels 
were  also  analyzed  for  each  feature.  Of  the  416  projects  generated  for 
the  104  features  at  McClellan  AFB,  only  309  were  saved  for  further 
analysis.  The  remaining  were  discarded  as  duplicates.  The  generated 
project  data  were  reduced  into  a  condensed  data  set  for  input  in  the 
Toyoda  algorithm.  The  condensing  program  took  less  than  4  central 
processing  seconds  to  extrapolate  the  important  data  from  the  dynamic 
programming  generated  listing  and  construct  a  data  matrix  for  the  Toyoda 
algorithm. 


Using  the  condensed  data  from  the  dynamic  program,  the  Toyoda 
program  optimally  allocates  a  pre-established  funding  limit  among  the 
candidate  projects  in  order  to  fund  those  projects  which  maximize  the 
TFW  weighted  structural  performance.  The  Toyoda  program  was  run  for 
fourteen  levels  or  network  funding  between  zero  and  25  million  present 
worth  dollars.  An  output,  similar  to  Figure  5.3  (Chapter  5),  was 
generated  for  each  funding  level.  Each  funding  level  output  and  the 
dynamic  programming  output  were  used  as  input  into  the  network  report 
program.  Examples  of  the  reports  are  shown  on  Figures  5.5,  5.6,  and  5.7 
(Chapter  5).  Using  the  information  produced  in  these  reports,  the 


consequences  of  varying  network  funding  level  can  be  made.  All  the 
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Figure  6.3  Saaple  of  Dynaoic  Progranning  Output  for  McClellan  AFB 


programs  developed  in  this  research  were  run  for  McClellan  AFB. 
Approximately  830  central  processing  seconds  were  used  for  the  base  as 
itemized  on  Table  6.1. 

6.3  McClellan  AFB  Application  Graph? 

Using  the  information  on  the  report  listings  from  the  various 
funding  levels,  the  consequences  of  specific  funding  levels  can  easily 
be  seen.  Figures  6.4,  6.5,  6.6,  and  6.7  are  constructed  from  data  in 
the  report  listings  at  each  level  of  network  funding  limit. 

Figure  6.4,  as  expected  has  a  smooth  continuous  increase  in  the 
objective  function  (TFW  weighted  structural  performance)  as  a  function 
of  network  funding  limit.  This  graph  amounts  to  relaxing  the  constraint 
of  spending  dollars  and  permits  more  projects  to  be  funded. 

Figure  6.5  plots  the  terminal  weighted  average  network  PCI  as  a 
function  of  network  funding  limit.  It  too,  has  a  continuous  increase  on 
terminal  PCI.  Note  that  the  first  time  the  curve  levels  off,  it  occurs 
at  a  PCI  corresponding  to  the  constraint  that  no  activities  will  be 
scheduled  until  the  pavement  reaches  that  level  (in  this  case,  70  PCI). 

Figure  6.6  plots  the  number  of  features  that  fall  below  a  PCI  of  40 
by  the  end  of  the  20  year  analysis  period.  There  is  a  continuous 
decrease  until  the  10  million  dollar  limit.  Analysis  of  the  data 


Table  6.1  Computer  Usage  Summary  for  McClellan  AFB  Application 


COMPUTER  USAGE  SUHMARY 
FOR  McCLEUAH  AFB 


Computer  Program  Central  Processing* 

Program  Description  Time  (Seconds) 


DYNC  Dynamic  Programing,  350 

generation  of  feature 
projects,  104  features 

CCNPRG  Condensing  Program,  4 

consolidates  OYNC  output 
for  Toyoda  Program 

TOY  Toyoda  Program,  375 

Network  allocation,  run 
14  times 

REPORT  Network  Report  Program  100 

compile  network  reports, 
run  14  times  _ 


N'ETWORK  FUNDS  SPENT,  MILLIONS,  PRESENT  WORTH 


139 


SUV3A  OS  JO  QN3  IV  lOd  1VNIWH31 


NETWORK  FUNDS  SPENT,  MILLIONS,  PRESENT  WORTH 


oinouioinoinotooino 
SHV3A  02  JO  <JN3  3H1  XV  0*  NVH1 


SS3T  lOd  HUM  S3HniV3J  JO  «30WnN 


NETWORK  FUNDS  SPENT,  MILLIONS,  PRESENT  WORTH 


141 


reveals  that  at  this  funding  level  all  features  with  any  value  to  the 
user  have  had  projects  funded.  With  additional  network  funding,  the 
Toyoda  algorithm  identifies  sore  gain  to  the  objective  function  by 
funding  more  expensive  projects  for  features  having  projects  selected  at 
lower  levels  of  funding. 

Figure  6.7  has  the  same  discontinuity  in  its  curve  as  it  plots 
network  performance  (animation  of  projects  performance)  versus  network 
funding  limit .  The  decrease  in  performance  with  increased  network 
funding  (after  10  million)  results  from  the  TFU  weighting  in  the 
objective  function. 

Figure  6.8  is  a  composite  summary  of  Figures  6.4,  6.5,  6.6  and  6.7 
shoving  on  one  graph  the  consequences  to  the  network  for  various  network 
funding  levels.  The  region  on  the  graph  between  5  and  10  million, 
indicates  a  network  funding  vhere  not  only  the  objective  function  levels 
off  but  the  other  three  network  consequence  indicators  also  level  off. 


4.4  Application  Discussion 


With  Jiese  graphs  and  reports,  airfield  pavement  engineers  can 
begin  to  effectively  manage  their  paveaent  system  (McClellan  AFB)  over 
the  analysis  period  (20  years)  for  major  rehabilitative  activities.  The 
graph  of  Figure  6.8  can  be  used  to  estimate  what  level  of  network 
funding  would  be  the  best  for  the  airfield  and  hence  they  can  determine 
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a  target  amount  to  request  free  the  funding  authorities. 


After  the  respasible  agency  finally  decides  on  the  network  funding 
limit,  the  actual  projects  to  be  funded  can  be  identified  frem  an 
existing  run  cf  the  Toyoda  algorithm  or  a  new  run  of  the  algorithm.  In 
either  case,  changes  to  decisions  made  at  the  network  level  can  be 
directly  reflected  to  specific  projects  that  comprise  the  network  plan. 

Just  as  important,  with  this  procedure  the  decision  makers  can  see 
the  consequences  of  their  decisions  in  quantifiable  terns  whether  their 
decision  is  to  reduce  or  increase  the  funding  limit  or  to  direct  that 
certain  activities  be  applied  to  specific  features. 


6.5  Summary 


This  chapter  discussed  the  application  of  the  procedures  developed 
in  this  research  to  an  actual  Air  Force  installation,  McClellan  AFB. 
Included  is  an  analysis  of  the  results  with  graphs  /plots  of  the  analysis 
data  to  illustrate  the  usefulness  of  the  procedures.  The  next  and  last 
chapter  includes  a  summary  with  r ecomaendations  and  conclusions. 
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CHAPTER  7 

SUMMARY,  RECOMMENDATION S ,  AMD  CONCLUSIONS 

7.1  General 

This  final  chapter  includes  a  stnsary  of  the  work  accomplished  in 
this  research,  recommendations  for  implementation  and  further  research, 
and  concluding  remarks  for  this  thesis. 

7.2  Sutrv  of  Research 

The  goal  of  this  research  has  been  to  develop  a  methodology  for 
managing  pavement  networks  over  prolonged  analysis  periods.  The 
developed  methodologies  are  flexible  and  responsive  to  decision  nakers 
and  engineers.  They  are  also  capable  of  reflecting  changes  in  the 
management  requirements  (reduced  budgets,  special  needs,  etc.)  on  not 
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only  a  macro  but  also  a  micro  scale.  In  other  words,  changes  to  the 
network  analysis  criteria,  are  reflected  in  specific  changes  in  which 
projects  are  selected. 

The  analysis  was  divided  into  two  parts,  project  level  and  network 
level  analysis.  Separate  methods  were  devised  for  each,  and  the  project 
level  procedures  were  designed  to  provide  inputs  into  the  network  level 
procedures.  But  each  level  of  analysts  can  still  be  used  independently 
of  each  other.  Althougn  the  project  level  procedures  are  tailored  to 
military  airfields,  the  network  analysis  procedures  are  general  enough 
to  be  used  with  most  resource  allocation  problms,  independent  of 
whether  or  not  they  are  pavement  related. 

The  project  level  analysis  procedures  optimally  select  and  schedule 
major  rehabilitation  activities  over  an  extended  analysis  period  (in 
this  research,  20  years).  The  activities  included  are:  routine 
maintenance,  reconstruction,  and  overlays  as  detailed  on  Chapter  3.  A 
computer  code  was  written  to  use  dynamic  programming  methods  to  select 
and  schedule  the  activities  over  the  analysis  period.  The  program 
optimized  the  activity  selection  and  scheduling  by  maximizing  the 
structural  performance  of  this  activity  for  the  pavement  section 
(feature).  The  structural  performance  was  characterized  by  the  area 
under  the  utility  weighted  PCI  (Pavement  Condition  Index)  versus  time 
(over  the  length  of  the  analysis  period)  curve.  The  information 
required  by  the  procedure  is  easily  obtained  at  each  Air  Base  and  no 
special  “pre-analysis'*  testing  or  investigation  is  necessary.  The 
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output  from  this  procedure  is  a  feature  by  feature  listing,  selecting 
and  scheduling  the  activities  that  maximize  the  structural  performance 
during  the  analysis  period  for  various  levels  of  funding  limits.  These 
results  can  be  used  by  decision  makers  and  pavement  engineers  to  help 
decide  vnich  rehabilitation  activities  should  be  considered  and  when  to 
schedule  them.  Hot  only  do  these  projects  have  a  present  worth  dollar 
cost,  but  they  also  are  given  a  quantified  measure  of  the  structural 
benefit  (performance)  of  the  project  with  respect  to  the  ±  rvement 
feature.  The  primary  purpose  of  developing  optimal  projects  for  each 
feature  was  to  generate  input  for  the  network  level  analysis. 

Once  a  set  of  projects  for  a  network  has  been  developed,  whether 
optimally  by  the  dynamic  programming  procedures  (recommended  way)  or  by 
other  methods  of  the  engineers'  choosing,  the  network  level  analysis 
method  can  be  used  to  select  the  best  of  these  projects  subject  to 
certain  constraints.  The  network  level  analysis  procedures  select,  from 
the  group  of  projects  submitted,  those  projects  that  maximized  the  sum 
of  the  user  value  weighted  structural  performance  of  each  project.  The 
mathematical  representation  of  this  selection  process  is  a  zero-one 
integer  linear  programming  model.  Although  solutions  to  a  large 
practical  size  problem  normally  require  an  is$ractical  aaount  of 
computation  time,  a  heuristic  developed  by  Toyoda  provides  extremely 
good  solutions. 

A  computer  code  was  written  to  transform  the  dynamic  programming 
generated  results  into  input  for  the  integer  programming  model,  and  then 
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Toycaa"s  method  is  used  to  select  those  projects  that  maximize  the  user 
weighted  (TFW)  structural  performance  (objective  function)  of  the  entire 
pavement  network.  Resulting  from  this  network  analysis  is  a  list  of 
selected  projects  (each  related  to  a  specific  feature)  that  maximizes 
the  objective  function  with  pre-established  constraints  (network  funding 
limit,  etc.)*  The  network  level  analysis  is  run  for  several  funding 
levels,  and  a  series  of  management  information  reports  are  generated  for 
each. 


With  these  reports,  the  consequences  of  selected  network  funding 
levels  can  quantitatively  be  compared.  In  addition,  an  estimate  of  an 
appropriate  level  of  funding  for  the  entire  system  can  be  made.  These 
reports  and  subsequent  analyses  also  reflect  the  consequences  of  network 
decisions  with  respect  to  specific  projects  affecting  individual 
features  in  the  pavement  network.  For  example,  the  impact  of  a  25Z 
reduction  in  network  funding  can  easily  be  seen  on  the  list  of  selected 
projects. 

The  entire  set  of  developed  procedures  and  implementing  programs 
were  developed  to  be  amenable  to  further  research  and  are  adaptable  to 
future  changes  (i.e.,  new  PCI  predictive  equations).  By  having 
implmented  them  within  modular  computer  programs,  these  procedures  can 
easily  incorporate  future  developments  in  pavement  technology. 

The  information,  data,  and  programs  for  this  thesis  research  are. 
available  through  Professor  Michael  I.  Darter  at  111  Talbot  Laboratory, 
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104  South  Wright  Street,  Orbana,  IL,  61801. 


7.3  Recommendations 


During  the  course  o£  this  research  several  questions  arose,  most 
were  handled  adequately  hut  others  vere  solved  in  an  expeditious  manner. 
Specifically,  those  areas  to  which  additional  research  could  be 
addressed  are: 

1)  In  this  work,  the  network  level  optimization  objective 
function  was  defined  to  be  the  structural  performance  of  the 
projects  to  the  respective  features  weighted  by  the  relative 
worth  of  the  projects'  affected  feature  to  the  user  (TFV) . 
This  is  calculated  as  the  siaple  product  of  TFW  of  the  feature 
and  the  total  area  under  the  utility  weighted  PCI  versus  time 
plot  for  the  project.  Use  of  this  parameter  for  optimization 
needs  further  research  to  validate  its  adequacy.  The  present 
method  is  a  compromise  between  an  objective  function  based 
solely  on  user  worth  of  the  feature,  or  solely  on  structural 
performance  of  the  project.  Neither,  by  themselves,  is 
correct,  but  vhich  formula  should  be  used  to  represent  the 
actual  combined  relationship  is  uncertain. 

2)  The  available  activities  to  the  dynamic  programming  algorithm 
should  be  increased  from  the  present  five  (routine  maintenance; 
reconstruction  of  PCC;  and  overlays  of  2,  4,  and  6  inches  of 
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AC)  to  include  such  activities  as:  recycling,  reconstruction 
with  flexible  pavement,  undersealing  and  grinding,  etc.*  Any 
of  these  activities  can  be  included  so  long  as  their  effects  to 
the  PCI  predictive  equation  can  be  quantitatively  measured. 

3)  These  programs  and  procedures  along  with  the  internal 
algorithms  need  to  be  validated  and/or  calibrated  to  fit  the 
decision  practices  of  current  Air  Force  engineers  and  managers. 
Where  opinion  and  judgement  were  necessary,  those  of  the  author 
were  used. 


7.4  Conclusions 


This  thesis  documents  research  conducted  over  the  past  28  months. 
During  the  course  of  this  work  both  project  level  and  a  network  level 
analysis  methodologies  have  been  developed.  This  research  is  the  first 
instance  in  which  that  operations  research  techniques  (dynamic  and 
integer  programming  algorithms),  has  been  applied  to  optimizing  the 
development  of  projects  for  individual  airfield  pavements  and  the 
selection  of  airfield  pavement  projects  for  an  entire  network  of 
pavement  sections. 

A  simple  comparison  ot  applying  traditional  manual  techniques  and 
those  techniques  developed  in  this  research  was  made,  and  an  application 
of  the  methods  to  a  real  Air  Force  base  wa6  presented.  The  application 


to  an  existing  Air  T  *ce  base  was  dene  to  show  the  feasibility  of 
applying  these  techniques.  What  renains  to  be  done  is  field  testing 
these  methods  to  incorporate  cuanges  and  further  validate  the 
procedures.  A  more  detailed  comparison  of  the  new  methods  and  the 
traditional  methods  of  pavement  engineers  would  be  difficult  because  the 
new  methods  provide  additional  insight  into  managing  a  network  of 
pavement  sections  (i.e.,  relationship  of  user  value  and  structural 
performance,  implications  of  a  staged  decision  process,  etc.). 

In  summary,  this  research  has  developed  new  methods  for  pavement 
management  which  are  implemen table  and  which  provide  better  pavement 
management  than  can  be  achieved  through  traditional  methods.  Future 
development  of  pavement  management  systems  for  both  hi2huavs  and 
airfields  must  begin  to  incorporate  methods  such  as  these  and  other  OR 
teenniques  in  order  to  assure  maximal  pavement  performance  under 
scarcities  of  resources  and  money. 
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